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ABSTRACT 


iii 


Extrinsic  Fabry-Perot  interferometric  (EFPI)  fiber  optic  sensors  were  used  to 
determine  the  first  five  modal  frequencies  of  laminated  glass/epoxy  composite  beams. 
EFPI  fiber  optie  sensors  and  piezoelectric  (PZT)  ceramic  sensor  both  found  identical 
modal  frequencies  of  the  composite  beams,  however,  EFPI  fiber  optic  sensors  showed 
more  sensitivity  and  better  signal-to-noise  ratios.  The  analytical  classical  beam  theory 
and  a  finite  element  model  validated  the  EFPI  modal  frequency  measurements.  Several 
8-ply  glass/epoxy  composite  beams,  each  26.04  cm  long  and  2.33  cm  wide  were 
fabricated.  Five  damaged  beams  with  specified  delaminations  as  well  as  two  different 
undamaged  beams  were  tested  for  modal  frequencies.  Five  different  delaminations  with 
sizes  ranging  from  1.27  cm  to  6.35  cm  long  were  incorporated  in  the  midplane  of  the 
beams.  Each  modal  frequency  shifted  with  a  change  of  delamination  size  and  location. 
Classical  beam  theory  was  used  to  simulate  modal  frequency  data  sets  for  1097  different 
prescribed  delamination  sizes  and  locations.  These  data  sets  were  applied  for  training  and 
testing  a  feedforward  backpropagation  neural  network.  Finally  the  neural  network  was 
tested  against  the  EFPI  fiber  optic  sensor  modal  frequency  results.  The  delamination  size 
and  location  predictions  resulted  from  neural  network’s  simulation  had  an  average  error 


of  5.9%  and  4.7%  respectively. 
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I.  INTRODUCTION 


Fiber  reinforced  composite  materials  are  being  used  more  frequently  in  many 
different  engineering  fields.  The  automobile,  aerospace,  naval,  and  civil  industries  all  use 
composite  materials  in  some  way.  Composite  materials  are  gaining  popularity  because  of 
high  strength,  low  weight,  resistance  to  corrosion,  impact  resistance,  and  high  fatigue 
strength.  Other  advantages  include  ease  of  fabrication,  flexibility  in  design,  and  variable 
material  properties  to  meet  almost  any  application  [1]. 

With  this  increase  in  structural  applications  for  composite  materials  comes  an 
increase  in  safety  awareness.  Unlike  steel  or  cement,  it  is  difficult  to  see  damage  to 
composite  materials  with  the  naked  eye.  The  delamination  of  laminated  composite 
structures  is  an  important  failure  mode  and  consequently  has  received  much  attention 
recently  [2].  Delaminations  affect  the  strength  and  integrity  of  the  composite  structure. 

Embedded  or  surface  mounted  sensors  can  be  used  to  continuously  monitor  the 
composite  material.  Real  time  measurements  of  strain,  temperature,  or  composite 
material  failure  can  be  determined  by  many  different  types  of  sensors;  fiber  optic, 
piezoelectric,  or  LVDT.  Long  term  health  monitoring,  fatigue  deformation,  and 
delamination  formation  can  also  be  monitored  with  sensors  [3]. 

Fiber  optic  sensors  hold  many  advantages  over  conventional  electrical  sensors. 
Fiber  optic  sensors  are  light  weight,  small  in  size,  have  high  sensitivity,  and  have  a  large 
bandwidth.  Fiber  optic  sensors  are  all-passive  (dielectric)  which  provides  immunity  to 
electromagnetic  interference  and  reduces  noise  [4].  Also,  small  size,  environmental 
ruggedness,  and  high  temperature  resistance  allow  fiber  optic  sensors  to  be  embedded 
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inside  the  composite  structure.  Conventional  electrical  sensors  cannot  be  imbedded  into 
the  composites  because  high  temperatures  change  their  performance  permanently  [2], 
Electrical  sensors  are  also  large  compared  to  fiber  optic  sensors  and  have  soldered 
electrical  connections. 

Fiber  optic  sensors  can  potentially  perform  four  functions  in  smart  structures 
applications.  First,  they  can  monitor  strain,  temperature,  and  pressure  during  the 
manufacturing  process.  Second,  they  perform  a  non-destructive  evaluation  of  structural 
integrity  at  any  point  during  the  manufacturing  process.  The  fiber  optic  sensors  can  also 
be  used  as  a  health  monitoring  system  for  the  composite  structure.  Finally,  they  can 
complement  performance  monitoring  and  control  systems  [3]. 

Smart  structures  is  a  term  used  to  describe  a  system  that  uses  a  sensor  or  network 
of  sensors,  a  processor,  and  actuators  to  sense  the  environment,  interpret  data,  make 
decisions,  and  respond  to  the  decision.  In  this  study  a  single  fiber  optic  sensor  and  a 
single  piezoelectric  lead  zironate  titanate  (PZT)  actuator  are  used  to  perform  real  time 
health  monitoring  of  a  glass/epoxy  composite  beam.  The  sensor  chosen  for  this 
experiment  is  an  extrinsic  Fabry-Perot  interferometric  (EFPI)  fiber  optic  sensor. 

In  this  research,  the  first  five  modal  frequencies  of  several  glass/epoxy  composite 
beams  with  and  without  various  prescribed  delaminations  are  measured  with  EFPI  fiber 
optic  sensors.  Classieal  beam  theory,  a  finite  element  model,  and  a  piezoelectric  (PZT) 
sensor  are  all  used  to  validate  the  accuracy  of  the  modal  frequencies  measured  using  the 
EFPI  fiber  optic  sensors.  The  classical  beam  theory  is  also  used  to  calculate  modal 
frequencies  for  variations  on  prescribed  delamination  size  and  location.  These  modal 
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frequencies  are  used  to  train  a  backpropagation  neural  network  to  predict  the  size  and 
location  of  the  prescribed  delaminations. 

The  next  section  is  a  review  of  the  literature  that  gives  background  information  as 
well  as  alternate  sources  for  the  theories  and  applications  discussed  in  this  study.  The 
third  section  is  an  overview  of  the  theoretical  concepts  used  to  measure  and  compare  the 
experimental  results.  Then  the  exact  procedures  used  to  complete  the  experiments  are 
outlined.  The  fifth  section  gives  the  results  of  the  experiments  followed  by  the  discussion 
of  these  results  in  the  sixth  section.  The  last  section  is  a  conclusion  of  this  study. 
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II.  LITERATURE  REVIEW 

A.  SMART  STRUCTURES  AND  SENSING 

A  “smart”  structure  is  one  which  incorporates  a  system  that  uses  a  sensor  or 
network  of  sensors,  a  processor,  and  actuators  to  sense  the  environment,  interpret  data, 
make  decisions,  and  respond  to  the  decisions.  The  development  of  a  full  smart  structures 
system  requires  the  coordination  of  many  different  engineering  disciplines.  The 
materials,  structures,  actuators,  sensors,  signal  processing,  and  systems  must  all  work 
together  to  make  smart  structures  a  reality. 

The  level  of  activity  in  smart  structures  research  has  greatly  increased  in  the  last 
decade  [2].  One  particular  area  of  interest  is  applying  smart  structures  to  monitor  real 
time  damage  assessment.  An  early  study  involved  simple  embedded  fiber  optic  grids  that 
located  damage  by  monitoring  loss  of  signal  from  the  optical  fibers  near  the  damage 
location  [5].  A  similar  experiment  was  conducted  on  a  composite  lattice  [6].  A  simple 
type  of  fiber  optic  sensor  involved  the  strain  induced  modulation  of  the  coherent 
interference  pattern  leaving  the  fiber.  This  was  the  basis  for  early  fiber  optic  technology 
[7].  Now  more  sensitive  techniques  using  optical  interference  have  been  developed. 
These  include  the  Fabry-Perot  interferometer,  the  Mach-Zehner  interferometer,  and  the 
Michelson  interferometer  [8].  The  suitability  of  fiber  optic  sensors  in  smart  structures 
applications  was  demonstrated  in  a  study  incorporating  solid  state  actuators  and  neural 


processing  [9]. 
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B.  USE  OF  HBER  OPTIC  SENSORS 

Fiber  optic  sensors  are  ideal  for  smart  structures  applications  and  for  use  with 
composite  materials.  They  have  low  weight,  small  size,  high  sensitivity,  and  a  large 
bandwidth.  Fiber  optic  sensors  are  all-passive  (dielectric)  which  provides  immunity  to 
electromagnetic  interference  and  reduces  noise.  Also,  small  size,  environmental 
ruggedness,  and  high  temperature  resistance  allow  fiber  optic  sensors  to  be  embedded 
inside  the  composite  structure. 

In  the  area  of  dynamic  testing,  the  extrinsic  Fabry-Perot  interferometric  (EFPI) 
fiber  optic  sensor  is  the  most  widely  studied.  A  Fabry-Perot  interferometer  as  a  fiber 
optic  sensor  was  first  introduced  in  1982  [10].  In  a  later  study  intrinsic  Fabry-Perot 
interferometric  fiber  optic  sensors  were  embedded  into  composite  material  and  used  as 
ultrasound  sensors  [11].  Next,  EFPI  fiber  optic  sensors  were  used  as  strain  sensors  for 
composite  and  metal  beams  as  well  as  for  composite  plates  [12,13,14,15,16].  Many 
signal  processing  techniques  and  sensor  configurations  were  also  discussed  in  these 
studies.  Fatigue  testing  was  done  on  a  full-scale  F-15  aircraft  and  the  EFPI  fiber  optic 
sensor  was  applied  to  smart  structure  applications  and  nondestructive  health  monitoring 
[17,18,19].  Modal  frequencies  and  mode  shapes  were  detected  using  fiber  optic  sensors 
in  1994  [20].  In  1995,  EFPI  fiber  optic  sensors  were  used  to  detect  delaminations  in 
composite  beams  [21].  References  20  and  21  both  relate  directly  to  the  present  study.  As 
an  alternative  to  EFPI  fiber  optic  sensors,  fiber  Bragg  gratings  and  birefringent  integrated 
fiber  optic  sensors  have  been  used  for  vibration  sensing  of  composite  structures  [22,23]. 
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C.  COMPOSITE  MATERIAL  ANALYSTS 

Delamination  detection  in  smart  composite  structures  has  received  much  attention 
recently  because  delaminations  are  an  important  failure  mode  for  laminated  composites. 
Delaminations  affect  the  strength  and  integrity  of  the  composite  structure  and  may  cause 
structural  failure  at  a  load  lower  than  the  design  load.  The  effects  of  delamination  on 
buckling  load,  postbuckling  deformation,  and  delamination  growth  under  various 
geometrical  parameters,  dynamic  and  static  loading  conditions,  material  properties,  and 
boundary  conditions  have  been  studied  extensively  for  the  past  few  years  [24,25].  These 
studies  predicted  the  size  and  location  of  delaminations  in  composite  structures. 
Interlaminar  fracture  toughness  based  on  static  models  have  been  used  to  predict  the  size 
of  delamination  cracks  in  composite  laminates  [26].  The  use  of  dynamic  testing  to  detect 
and  quantify  the  effects  of  internal  damage  on  modal  parameters  has  been  studied  on  a 
limited  basis  over  the  past  twenty  years. 

Shifts  in  natural  frequencies  have  been  used  to  identify  damage  generated  by 
various  mechanical  loadings  including  prescribed  delamination  [27].  Effects  of 
prescribed  delamination  position  and  size  on  the  natural  frequencies  of  the  graphite  epoxy 
beam  specimens  have  been  investigated  analytically  and  experimentally  [28,29].  The 
analytical  approach  was  based  on  classical  beam  theory  and  imposed  the  constraint  of 
identical  transverse  displacements  for  the  upper  and  lower  sublaminates.  A  finite  element 
model  has  been  utilized  to  detect  and  suppress  delamination  based  on  such  piezoelectric 
behaviors  [30].  Piezoelectric  sensors  have  been  used  to  measure  the  natural  frequencies  of 
composite  beams  before  and  after  prescribed  delaminations  and  therefore  indicate  the 
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presence  and  size  of  the  delamination  [31],  These  results  were  compared  to  theoretical 
predictions  of  an  in-house  finite  element  model  based  on  shear  deformation  theory  [32], 
Fiber  optic  sensors  have  been  used  to  measure  acoustic  signatures,  changes  in 
strain  profiles,  delamination,  and  other  indications  of  changes  in  structural  characteristics 
of  fabricated  composite  structures  [6,11,16,21].  Fiber  optic  sensors  have  found 
increasing  applications  in  development  of  smart  composite  structures  because  of  their 
superior  advantages  including  compatibility  with  their  host  structures  and  the  ability  to  be 
easily  multiplexed  [9]. 


D.  NEURAL  NETWORK  APPLICATIONS 

Neural  networks  have  been  used  as  effective  tools  to  predict  and  generalize 
unknown  parameters  in  physical  systems  [33].  Recently  neural  network  simulations  have 
been  accurate,  robust,  and  fast  in  response  to  changes  in  dynamic  characteristics  of 
composite  structures.  In  particular,  they  have  been  utilized  in  such  applications  as 
damage  assessment  and  fatigue  monitoring  of  composite  structures  [34,35].  Neural 
networks  have  joined  feedback  control  circuits  and  programmable  logic  controllers  as 
processing  elements  in  smart  structures.  Neural  networks  are  especially  suited  for 
applications  in  which  large  numbers  of  sensors  and  actuators  are  used.  Neural  networks 
have  parallel  computing  architecture  that  allows  multiple  inputs  and  outputs  to  be 
processed  in  nanoseconds.  Also,  neural  networks  can  learn  from  experimental  data  and 
adapt  to  changing  conditions.  This  fast  processing  speed  and  adaptability  provide  many 
possible  applications  for  neural  networks  in  smart  structures. 
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III.  THEORY 

A.  SENSORS  AND  ACTUATORS 

Extrinsic  Fabry-Perot  interferometric  (EFPI)  fiber  optic  sensors  are  the  primary 
sensors  used  in  this  study.  EFPI  fiber  optic  sensors  are  used  to  experimentally  measure 
the  modal  frequencies  of  a  cantilever  glass/epoxy  composite  beam.  Piezoelectric  devices 
are  also  used  as  sensors  and  actuators.  A  piezoelectric  (PZT)  ceramic  sensor  is  used  on 
three  different  beams  to  verify  the  accuracy  of  the  EFPI  fiber  optic  sensors.  A  PZT 
ceramic  actuator  is  used  on  all  seven  beams.  The  PZT  ceramic  actuator  is  used  to  vibrate 
the  beam  at  different  frequencies  so  a  frequency  response  can  be  measured  for  the  beam. 

1.  Fabry-Perot  Interferometer.  The  Fabry-Perot  interferometer  consists  of  two 
mirrors  of  reflectance  R,  and  Rj  separated  by  a  distance  L.  The  reflectance  of  the  system 
is  given  by  the  ratio  of  reflected  power  P^  to  the  incident  power  Pj.  Assuming  the  mirrors 
are  of  equal  reflectivity  (R^  =  R^  =  R): 

_  2R(l-cos0) 

l  +  R-2R^(cos0)  ^  ’ 

where  ([)  is  the  round  trip  phase  shift  of  the  coherent  light  inside  the  cavity  [2].  The  phase 
shift  is  given  by: 

(t)  =  47rnL/^  (2) 

where  L  is  the  cavity  length,  n  is  the  index  of  refraction  of  the  cavity,  and  X  is  the  free 
space  wavelength  of  the  coherent  light  source.  For  this  study,  the  cavity  media  is  air  (n  = 
1.0)  and  the  laser  diode  wavelength  is  1310  nm.  The  finesse  F  is  a  parameter  used  to 
describe  the  sensitivity  of  the  interferometer.  Finesse  is  defined  as  the  ratio  of  the  phase 
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change  between  adjacent  transmittance  peaks  to  the  phase  change  between  half-maximum 
points  on  either  side  of  a  peak  [4].  It  can  be  written  as: 


F  = 


AR 

{\-Rf 


(3) 


Therefore,  as  R  increases,  so  does  F. 

2.  Extrinsic  Fabrv-Perot  Interferometric  Sensors.  An  EFPI  fiber  optics  sensor  is 
schematically  shown  in  Figure  1(a).  The  EFPI  fiber  optic  sensor  utilizes  multiple-beam 
interference  between  two  polished  end  faces  of  a  single  mode  fiber  and  a  multimode 
fiber.  The  Fresnel  reflection  from  the  glass-air  interface  at  the  front  of  the  air  gap 
(reference  reflection)  and  the  reflection  from  the  air-glass  interface  at  the  far  end  of  the  air 
gap  (sensing  reflection)  interfere  in  the  single  mode  fiber.  The  far  end  of  the  multimode 
fiber  is  shattered  so  the  reflections  from  the  far  end  do  not  add  to  the  detector’s  noise. 
The  phase  of  the  light  between  the  two  endfaces  in  the  air  gap  region  is  modulated  in 
response  to  vibration  via  the  strain-optic  effect.  The  active  sensing  length,  or  gauge 
length,  is  length  of  the  capillary  tube  that  is  cemented  around  the  two  fibers. 

The  EFPI  fiber  optic  sensors  used  in  this  study  had  gauge  lengths  of  8mm  and 
cavity  lengths,  L,  of  100  pm.  They  also  had  low  finesse,  (R  «  1)  so  equation  (1)  can  be 
approximated  by; 

^=2R(l-cos<t)).  (4) 


The  sensor  detects  vibration  by  changes  in  (j)  as  a  function  of  L  (see  equation  2).  As 
vibration  modulates  the  axial  sensor,  strain  changes  proportionally.  This  change  in 
reflected  power  is  detected  by  the  fiber  optic  support  system. 
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The  reflected  optical  fringe  response  of  the  EFPI  fiber  optic  modal  frequency  sensor  is 
converted  into  electrical  signals  by  a  fast  response  detector.  The  modal  frequency  content 
of  the  electrical  signal  is  obtained  by  a  Fourier  transform  unit,  which  is  a  spectrum 
analyzer  is  used  in  this  study. 

3.  Piezoelectric  Ceramic  Sensors  and  Actuators.  The  piezoelectric  effect  refers  to 
the  conversion  of  applied  mechanical  stress  in  a  piezoelectric  material  into  electric  field 
[3].  An  applied  electric  field  can  also  produce  mechanical  stress  in  the  material.  Sensing 
is  based  on  the  first  case  and  actuation  on  the  latter.  A  material  must  be  anisotropic  and 
electrically  poled  in  order  for  the  piezoelectric  effect  to  occur  [3]. 

The  response  of  a  piezoelectric  material  is  determined  by: 

D  =  dX  +e^e„E  (5) 

and 


X  =  sX  +dE  (6) 

where  D  is  the  electric  displacement  vector,  E  is  the  electric  field  vector,  is  the  relative 
dielectric  constant,  8„  is  the  free  space  dielectric  constant,  X  is  the  stress  applied  to  the 
material,  x  is  the  strain,  s  is  the  mechanical  compliance  at  a  particular  value  of  E,  and  d  is 
the  piezoelectric  strain  coefficient.  D,  X,  x,  and  E  are  vectors  and  d,  s,  are  tensors.  By 
rearranging  equations  (5)  and  (6)  and  assuming  D  is  constant,  the  strain  becomes: 

X={d-[(s8,£„)/d]}E.  (7) 

In  this  study,  identical  PZT  ceramic  patches  were  used  as  sensors  and  actuators. 
However,  a  major  difference  exists  between  using  PZT  ceramic  material  as  a  sensor  and 
as  a  actuator.  PZT  ceramic  material  converts  strain  into  small  electric  fields,  typically 
kilovolts/meter,  for  sensor  applications.  For  actuator  applications  the  electric  field  must 
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be  very  large,  approximately  1  megavolt/meter,  to  produce  a  mechanical  strain  (50  GPa 
to  100  GPa).  Typically  strain  levels  that  can  be  applied  by  PZT  ceramic  actuators  are 
0.1%. 

B.  COMPOSITE  MATERIAL  FABRIC ATTON 

Two  27.94  cm  x  27.94  cm  glass/epoxy  composite  plates  were  fabricated  in  the 
Composite  Fabrication  Laboratory  at  the  University  of  Missouri-Rolla.  The  plates  were 
eight  layer  symmetric  laminates  [0/90/0/90]5.  The  lamina  orientation  can  be  seen  in 
Figure  2.  The  material  properties  of  the  composite  plates  are  the  following:  £„  =  42.34 
GPa,  Ejj  =  1 1.72  GPa,  Gjj  =  7.10  GPa,  v,2  =  0.27,  and  p  =  1901.5  kg/m^  where  Ej,  and  E^^ 
are  the  principal  moduli,  G^^  is  the  shear  modulus,  v,^  is  the  Poisson  ratio,  and  p  is  the 
density.  The  plates  were  fabricated  using  a  75  ton  Drake  Model  #44-806  hydraulic  hot 
press.  The  fabrication  curve  used  to  control  the  hot  press  during  plate  fabrication  can  be 
seen  in  Figure  3.  The  plain  solid  line  is  the  manufacturer’s  recommended  fabrication 
curve.  The  squares  are  the  actual  temperatures  measured  by  the  hot  press  during  plate 
fabrication.  A  diamond  saw  was  used  to  cut  the  first  plate  into  beams  27.94  cm  long, 
2.54  cm  wide,  and  0.18  cm  thick.  Teflon  strips,  0.005  cm  thick,  were  inserted  in  the 
midplane  of  the  second  composite  plate  before  fabrication.  The  delamination  placement 
can  be  seen  in  figure  4.  The  plate  was  then  cut  into  ten  beams  27.94  cm  long,  2.33  cm 
wide,  and  0.15  cm  thick.  The  beams  contained  0  cm,  1.27  cm,  2.54  cm,  3.81  cm,  5.08 
cm,  and  6.35  cm  delaminations,  respectively,  located  at  the  center  of  the  beam. 


Figure  3.  Heating  and  cooling  curve  for  the  glass/epoxy  composite  plate  fabrication. 


C.  COMPOSITE  MATERIAL  ANALYvSTS 

Classical  beam  theory  approximation  is  a  one-dimensional  model  used  to  evaluate 
the  effects  of  delamination  position  and  size  on  the  beam  natural  frequencies  [28].  The 
formulation  is  based  on  a  constrained  vibration  model  with  the  assumption  of 
infinitesimal  and  identical  transverse  displacements  of  the  upper  and  lower  sublaminates. 
The  model  allows  for  independent  extensional  and  bending  stiffness  and  does  not  include 
bending/extensional  coupling  stiffness  terms.  Therefore  this  model  works  best  for 
symmetric  laminates.  Figure  4  represents  the  delamination  model  of  a  beam  of  unit 
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width,  thickness  h,  and  length  L.  Delamination  is  assumed  to  be  uniform  across  the 
width  of  the  beam.  The  delamination  has  length  and  goes  from  x  =  Ij  to  x  =  Ij  +  1^ 
parallel  to  the  neutral  axis  and  at  a  distance  of  h^  from  the  upper  surface  of  the  beam. 


Figure  4.  Cantilever  composite  beam  with  prescribed  delamination  location.  Segments 
land  4  are  the  base  laminates  and  segments  2  and  3  are  sublaminates  above  and 
below  the  prescribed  delamination. 


The  delaminated  beam  is  divided  into  four  sections  (see  figure  4)  and  each 
segment  is  characterized  by  an  extensional  stiffness  A,  a  bending  stiffness  D,  and  a  mass 
density  per  unit  length  m.  The  governing  differential  equations  for  eaeh  segment  are  the 
following: 


D.  4  -  m.(0  w  =  0 
dx 


(5) 


and 

dP.  d^U; 

— -  =  A  - -  =  n 

dx  '  dx^ 


(6) 
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where  Uj  and  Wj  are  axial  and  transverse  deflections  and  Pj  is  the  axial  force  in  each 
segment  (note  Pj  =  =  0).  The  solution  to  equations  (5)  and  (6)  has  the  form: 

Wj  =  C|COsh(kx)  +  Cj^sinhCkx)  +  Cj3Cos(kx)  +  Cj4sin(kx)  (7) 

and 

Ui  =  Bio  +  B„x  (8) 

where  k  is 


m-CD^ 


(9) 


for  i  =  1,2, 3, 4.  There  are  20  unknown  constants  in  Equations  (7)  and  (8).  However,  the 
deflections  of  segments  2  and  3  are  forced  to  vibrate  together  by  setting  C^,  =  Cjj,  = 
C33,  C23  =  C33,  and  €34  =  C34.  These  constraints  are  for  analysis  of  midplane  delaminations. 
Therefore  only  16  boundary  and  continuity  conditions  are  needed  for  the  16  remaining 
unknowns  in  Equations  (7)  and  (8).  For  a  cantilever  beam,  one  end  of  the  beam  is 
clamped  and  the  other  end  is  freely  vibrating.  The  16  geometric  boundary,  continuity, 
moment  balance,  shear  force  balance,  and  compatibility  for  axial  midplane  displacements 
conditions  are  thereby  set.  These  conditions  yield  16  homogeneous  equations  with  16 
unknowns.  The  16  coupled  equations  can  be  seen  in  Appendix  A,  program  Delam.m.  In 
matrix  form  the  equations  are: 

[Q,(k)][V.]  =  [0]  (10) 

where  Qij(k)  are  the  coefficients  of  the  unknown  constants  and  V.  are  the  elements  of  a 
column  vector  containing  the  unknown  constants.  The  expressions  for  the  256  constants 
in  Equation  (10)  are  given  elsewhere  [29].  For  a  non-trivial  solution,  the  determinant  of 
the  coefficient  matrix  Qyfk)  must  equal  zero.  Using  MATLAB  the  expressions  are 
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symbolically  derived  for  the  256  constants  as  a  function  of  k  and  other  beam  parameters. 
Numerically  the  determinant  of  Qij(k)  is  calculated  for  different  values  of  k.  If  there  is  a 
sign  change  in  the  determinant  of  Qy(k)  for  two  adjacent  values  of  k,  then  the  determinant 
of  Qi|(k)  has  a  root  in  between  these  two  adjacent  values.  The  root  is  accurately 
determined  by  repeated  bisection  of  the  interval  between  the  two  adjacent  k  values 
(Appendix  A,  program  Final.m).  The  modal  frequencies  are  then  found  using  equation 
(9).  The  bending/extensional  coupling  is  not  allowed  in  this  model. 

The  finite  element  model  discussed  in  this  study  was  used  to  determine  the 
theoretical  modal  frequencies  for  an  undamaged  glass/epoxy  composite  beam.  Free 
vibrations  of  clamped  composite  beams  based  on  the  classical  theory  which  neglects  the 
effects  of  transverse  shear  deformation  and  rotary  inertia  have  been  considered  to 
determine  the  five  modal  frequencies  of  a  glass/epoxy  beam  [36,37].  For  a  symmetric 
cantilever  composite  beam  the  equation  for  the  modal  frequencies  is: 


/  =  — 
J  n 


KL  R 


2kL  \  pA 


(11) 


where  L  is  the  length,  A  is  the  cross-sectional  area,  p  is  the  mass  density,  and  Dj,  is  the 
flexural  stiffness  of  the  laminated  composite  beam.  The  mode  parameter  k„  satisfies  the 
following  equation: 

cos(k„L)  cos(k„L)  +1=0.  (12) 

First-order  shear  deformation  and  rotary  inertia  should  be  included  for  more 


accurate  vibration  analysis  of  laminated  composite  beams  [38].  In  this  study  the  results 
from  a  finite  element  model  based  on  a  higher  order  shear  deformation  theory  are  also 
used  [32].  Results  can  be  seen  in  Appendix  C.  The  finite  element  model  accounts  for  the 
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effects  of  in-plane  inertia  and  rotary  inertia  in  the  formulation  of  the  mass  matrix.  Also 
the  Poisson  effect  is  incorporated  in  the  formulation  of  the  beam  constitutive  equation. 

Lagrange  linear  interpolation  functions  and  Hermite  cubic  interpolation 
functions  are  used  to  interpolate  the  generalized  displacements.  The  element  equation 
after  assembling  the  elements  becomes: 

([K]-(D^[M]){A}=0  (13) 

where  [K]  is  the  element  stiffness  matrix,  [M]  is  the  mass  matrix,  ro  is  the  angular 
frequency,  and  {A}  is  the  mode  shape.  Equation  (13)  can  be  solved  after  imposing  the 
boundary  conditions.  The  natural  frequencies  of  a  cantilever  composite  beam  are  the 
object  in  this  study.  Therefore  the  boundary  conditions  are; 

u  =  0,  w  =  0,  (t)  =  0  (14) 

for  the  clamped  edge.  Using  the  results  from  a  Fortran  finite  element  code,  natural 
frequency  solutions  are  obtained  for  the  free  vibration  of  laminated  composite  beams. 
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IV.  EXPERIMENTAL  PROCEDURES 

A.  FREQUENCY  RESPONSE  TESTING 

Two  27.94  cm  x  27.94  cm  glass/epoxy  composite  plates  were  fabricated  with  an 
eight  layer  symmetric  laminate  structure  [0/90/0/90]jg.  The  first  plate  was  cut  into  beams 
27.94  cm  x  2.54  cm  x  0.18  cm  which  contained  no  delaminations.  Teflon  strips,  0.005 
cm  thick,  were  inserted  in  the  midplane  of  the  second  composite  plate.  This  plate  was 
then  cut  into  ten  beams  27.94  cm  long,  2.33  cm  wide,  and  0.15  cm  thick.  The  beams 
contained  0  cm,  1.27  cm,  2.54  cm,  3.81  cm,  5.08  cm,  and  6.35  cm  delaminations 
respectively  located  at  the  center  of  the  beam. 

Figure  5  shows  the  experimental  setup.  A  PZT  ceramic  patch,  for  use  as  an 
actuator,  was  bonded  to  the  underside  of  the  each  beam  3.81  cm  from  one  end.  Eccobond 
57  C  Conductive  Epoxy  was  used  to  bond  the  PZT  ceramic  patch  to  copper  foil.  The  foil 
was  then  attached  to  the  beam  using  general  purpose  Super  Glue.  A  second  PZT  ceramic 
patch,  for  use  as  a  sensor,  was  bonded  in  the  same  manner  2.54  cm  from  the  opposite  end 
of  two  different  undamaged  beams  and  the  5.08  cm  delaminated  beam.  An  extrinsic 
Fabry-Perot  fiber  optic  sensor  (EFPI)  was  mounted  2.54  cm  from  the  end  of  each  beam 
using  Super  Glue.  Each  beam  was  clamped  1.91  cm  from  the  end  closest  to  the  actuating 
PZT  ceramic  patch.  The  cantilever  beam  length  was  26.04  cm. 

An  HP35665A  spectrum  analyzer  was  programmed  to  produce  a  sine  sweep 
ranging  from  1  Hz  to  1000  Hz,  with  an  entry  step  size  of  0.1  Hz,  a  resolution  of  801 
points  per  sweep,  and  a  source  voltage  of  1  V^.  Trial  and  error  were  used  to  determine  the 
settings  that  produced  the  best  results.  The  output  signal  was  amplified  with  a  Piezo 


Systems,  Inc.  piezo  amplifier,  model  60.101.  The  amplifier  was  connected  to  the 
actuating  PZT  ceramic  patch.  The  EFPI  fiber  optic  sensor  was  connected  to  a  fiber  optic 
support  system  using  a  butt  coupler  and  index  matching  gel.  The  fiber  optic  support 
system,  model  FOSS- 1300-2,  was  manufactured  by  the  F&S  Fiber  Sensor  Technologies. 
The  fiber  optic  support  system  contains  a  laser  diode  (^=1300  nm),  an  optical  isolator,  a  3 
dB  fiber  coupler,  and  a  high  speed  photodetector  as  seen  in  Figure  1.  The  fiber  optic 
support  system  and  the  PZT  ceramic  sensor  output  signals  were  connected  to  the  inputs 
of  the  spectrum  analyzer  (see  Figure  5).  The  actual  experimental  setup  can  be  seen  in  the 
photograph  in  Figure  6. 


27.94  cm 


Figure  5.  Experimental  setup  for  determining  modal  frequencies  in  a  cantilever  composite 
beam. 
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Figure  6.  Actual  photograph  of  the  experimental  setup.  The  composite  beam  containing  a 
PZT  ceramic  actuator  and  a  EFPI  fiber  optic  sensor  is  at  the  far  left  of  the 
picture.  Next  to  the  beam  is  the  fiber  optic  support  system  (F&S  box).  The 
spectrum  analyzer  is  on  the  table  with  the  piezo  amplifier  on  top  of  it.  The 
plotter  is  to  the  far  right. 

The  sine  sweep  was  initiated  using  the  spectrum  analyzer  and  the  amplifier.  The 
output  signals  from  the  fiber  optic  support  system  and  the  PZT  ceramic  patch  were 
collected  by  the  spectmm  analyzer  in  the  frequency  domain,  and  the  first  five  modal 
frequencies  were  determined  and  plotted.  Both  undamaged  and  delaminated  glass/epoxy 
cantilever  beams  were  tested  and  the  modal  frequencies  were  successfully  found.  The 
experimental  modal  frequencies  were  verified  with  calculated  modal  frequencies  obtained 
from  classical  beam  theory  and  from  the  finite  element  model. 
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B.  NEURAL  NETWORK  SETUP 

A  backpropagation  neural  network  with  one  input  layer,  three  hidden  layers,  and 
one  output  layer  was  used  to  predict  delamination  size  and  delamination  location.  There 
are  no  established  procedures  for  choosing  the  optimal  number  of  hidden  layers  and 
processing  elements  per  layer.  Incorporation  of  more  processing  elements  enhances  the 
generalization  capability  of  the  network  for  a  larger  number  of  training  data  points.  The 
input  layer  has  five  processing  units  and  each  corresponds  to  one  of  the  first  five  modal 
frequencies  calculated  using  MATLAB  simulations  based  on  classical  beam  theory.  The 
first  hidden  layer  has  four  processing  elements  (neurons),  the  second  hidden  layer  has 
five  processing  elements,  and  the  third  hidden  layer  has  four  processing  elements.  The 
output  layer  consists  of  two  processing  elements.  The  network  predicts  delamination  size 
and  delamination  location.  NeuralWorks  Professional  II  software  was  used  to  design, 
train,  and  test  the  neural  network  on  a  Pentium  133  MHz  personal  computer.  A  basic 
schematic  of  the  neural  network  is  shown  in  Figure  7. 

In  this  study  the  hyperbolic  tangent  function  was  used  as  the  activation  (transfer) 
function.  Training  consisted  of  providing  a  set  of  1066  calculated  input  modal 
frequencies  and  output  delamination  size  and  delamination  location  to  the  network.  Each 
input  set  consisted  of  five  modal  frequencies  for  each  delamination  size.  The 
delamination  sizes  ranged  from  zero  to  20.22  cm.  There  were  two  to  ten  delamination 
locations  for  each  delamination  size.  The  number  of  delamination  locations  depended  on 
delamination  size  because  delaminations  cannot  exceed  the  beam  length.  The  network 
iteratively  adjusted  the  weights  of  each  of  the  nodes  to  obtain  the  desired  outputs.  The 
normalized  cumulative  delta  rule  is  used  to  train  the  network.  During  the  backward 
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propagation,  weights  in  nodes  were  adjusted  based  on  the  errors  using  the  gradient 
descent  equation  [39]. 


Input  layer  Hidden  layer  1  Hidden  layer  2  Hidden  layer  3  Output  layer 

Figure  7.  A  feed  forward,  backpropagation  neural  network  with  five  inputs  and  two 
outputs. 

To  optimize  the  performance  of  the  network,  learning  rate  and  momentum 
coefficients  can  be  adjusted.  Among  the  neural  network  parameters,  learning  rate  and 
momentum  coefficients  are  the  critical  variables  that  effect  the  convergence  of  the 
network.  The  learning  rate  started  at  0.5  and  was  reduced  as  iterations  increased  and  the 
momentum  coefficient  started  at  0.3.  The  initial  learning  rate  and  momentum  coefficients 
were  chosen  after  several  trials  to  find  optimum  performance. 

At  the  beginning  of  the  training,  the  network  produced  outputs  which  most  likely 
differ  from  the  desired  values.  The  procedure  was  repeated  until  the  difference  between 
the  desired  and  actual  outputs  was  reduced  to  a  preset  convergence  limit.  The  training 
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was  stopped  after  160,000  cycles.  The  root-mean-squared  error  converged  to  0.0005  at 
approximately  100,000  cycles.  After  100,000  cycles,  the  root-mean-squared  error  varied 
from  0.0005  to  0.1.  The  root-mean-squared  error  stayed  at  approximately  0.0005  until 
single  error  spikes  of  0. 1  would  appear.  The  spikes  were  nonperiodic  and  did  not  occur 
very  often.  These  spikes  were  assumed  to  be  point  errors  and  were  neglected.  Therefore, 
the  system  converged.  The  network  simulation  took  about  90  seconds  to  complete  on  a 
Pentium  133  MHz  personal  computer. 
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V.  EXPERIMENTAL  RESULTS 

A.  COMPSITE  BEAM  MODAL  FREQUENCIES 

Figure  8  shows  the  experimental  frequency  response  of  an  undamaged  composite 
beam.  The  response  functions  were  generated  using  a  curve  fitting  procedure  built  into 
the  HP35665A  spectrum  analyzer.  Each  peak  corresponds  to  a  modal  frequency  of  the 
beam.  Figure  8(a)  is  the  response  from  the  EFPI  fiber  optic  sensor  while  Figure  8(b)  is 
the  response  of  the  PZT  ceramic  sensor. 

As  described  earlier,  both  classical  beam  theory  and  a  finite  element  method  were 
used  to  calculate  modal  frequencies  of  the  undamaged  test  beam.  Table  I  compares  the 
experimental  frequencies  with  those  obtained  from  the  classical  theory  and  the  finite 
element  model.  Experimental  modal  frequencies  are  in  better  agreement  with  the 
calculated  modal  frequencies  from  the  classical  beam  theory  than  with  those  from  the 
finite  element  model.  The  calculated  difference  in  modal  frequencies  between  the  finite 
element  model  and  the  classical  beam  theory  is  also  shown  in  Table  I.  For  high  modes 
there  is  less  difference  between  the  modal  frequencies  calculated  from  the  finite  element 
model  and  classical  beam  theory. 

The  effect  of  the  beam  width  on  modal  frequencies  was  also  observed.  Two 
undelaminated  beams,  one  having  a  width  of  2.54  cm  and  the  other  having  a  width  of 
2.33  cm,  were  tested.  As  stated  earlier.  Figure  8  is  the  frequency  response  of  the  2.33  cm 
wide  beam.  Figure  9  is  the  frequency  response  for  the  2.54  cm  wide  beam.  Again,  the 
frequencies  are  identical  for  the  EFPI  fiber  optic  sensor  and  the  PZT  ceramic  sensor.  The 
peaks,  or  modal  frequencies,  of  the  two  different  beams  shifted  with  the  change  in  width. 


26 


Table  II  shows  the  result  of  the  modal  testing.  The  average  difference  between  the  modal 
frequencies  of  the  two  beams  is  12.9%.  However,  the  average  difference  in  the  finite 
element  code  between  a  2.54  cm  beam  and  a  2.33  cm  beam  is  26.2%.  Also,  the  average 
difference  in  the  modal  frequencies  of  the  two  beams  calculated  by  the  classical  beam 
theory  is  20.7%. 


Table  I:  Experimental  modal  frequencies  of  an  undamaged  beam  compared  to  MATLAB 
generated  frequencies  of  classical  beam  theory  and  frequencies  calculated  using  a 
finite  element  model  (FEM). 


Modal 

Frequencies 

Experiment 

(Hz) 

Classical 
Theory  (Hz) 

Difference 

Exp-Theory 

Finite 

Element  (Hz) 

Difference 

Exp-FEM 

Difference 

Theory-FEM 

1 

13.5 

16.7 

19.2% 

15.3 

11.7% 

8.4% 

2 

88.4 

100.8 

12.3% 

96.2 

8.1% 

4.5% 

3 

250.8 

282.1 

11.1% 

270.8 

7.4% 

4.0% 

4 

491.8 

552.9 

11.0% 

534.9 

8.1% 

3.3% 

5 

815.2 

913.9 

10.8% 

893.1 

8.7% 

2.3% 

Table  11:  Comparison  of  modal  frequencies  of  two  composite  beams  having  different 
widths. 


Modal 

Frequencies 

2.54  cm  Wide 
Beam  (Hz) 

2.33  cm  Wide 
Beam(Hz) 

Difference 

Difference 
in  FEM 

Difference 
in  Theory 

1 

12.2 

13.5 

9.6% 

26.1% 

23.3% 

2 

75.9 

88.4 

14.1% 

26.2% 

20.1% 

3 

213.3 

250.8 

14.9% 

26.2% 

20.1% 

4 

425.5 

491.8 

13.5% 

26.2% 

20.1% 

5 

714.0 

815.2 

12.4% 

26.2% 

20.1% 

Amplitude  (dB)  Amplitude  (dB) 


Figure  8.  Frequency  response  of  (a)  an  EFPI  fiber  optic  sensor  and  (b)  a  PZT  ceramic 
sensor  mounted  on  an  undamaged  glass/epoxy  beam.  The  beam  is  2.33  cm 
wide. 


Figure  9.  Frequency  response  of  (a)  an  EFPI  fiber  optic  sensor  and  (b)  a  Pzt  ceramic 
sensor  mounted  on  an  undamaged  glass/epoxy  beam.  The  beam  is  2.54  cm 
wide. 
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B.  DELAMINATED  COMPOSITE  BEAM  MODAL  FREQUENCIES 

The  five  delaminated  composite  beams  were  tested  after  the  two  undamaged 
beams.  Figure  10  shows  the  frequency  response  of  a  EFPI  fiber  optic  sensor  and  a  PZT 
ceramic  sensor  mounted  on  beam  containing  a  5.08  cm  delamination  incorporated  into 
the  midplane.  Again  the  modal  frequencies  were  identical  for  both  sensors.  Figure  1 1 
through  Figure  14  show  the  frequency  responses  of  beams  with  1.27  cm,  2.54  cm,  3.81 
cm,  and  6.35  cm  delaminations  respectively. 

The  experimental  values  for  the  modal  frequencies  are  found  in  Table  III. 
Examination  of  Figure  8  and  Figures  10  through  14  reveal  the  modal  frequencies  of  the 
delaminated  beams  are  shifted.  Table  III  summaries  these  shifts  in  modal  frequencies  for 
the  2.33  cm  wide  undamaged  beam  and  the  five  different  delaminated  beams. 


Table  III:  Experimental  modal  frequencies  of  delaminated  composite  beams. 


Modal 

Frequencies 

Undamaged 

(Hz) 

1.27  cm 
Delam  (Hz) 

2.54  cm 
Delam  (Hz) 

3.81  cm 
Delam  (Hz) 

5.08  cm 
Delam  (Hz) 

6.35  cm 
Delam  (Hz) 

1 

13.5 

17.2 

15.9 

15.9 

13.5 

14.7 

2 

88.4 

93.4 

95.9 

98.4 

94.7 

94.7 

3 

250.8 

256.9 

260.7 

254.5 

233.4 

210.8 

4 

491.8 

502.9 

515.5 

527.9 

519.2 

490.5 

5 

815.2 

835.1 

818.9 

750.3 

684.1 

621.6 

The  experimental  modal  frequencies  for  each  delamination  size  listed  in  Table  III 


were  compared  to  the  calculated  classical  beam  theory  modal  frequencies  found  in 
Appendix  B.  Table  IV  show  the  results  of  this  comparison.  A  percent  difference  was 
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also  calculated.  The  percent  difference  ranges  from  zero  to  16.5%  with  an  average 
difference  of  5.6%. 


Table  IV;  Comparison  of  experimental  and  classical  beam  theory  modal  frequencies  for 
different  delamination  sizes. 


Delamination 
Size  (cm) 


1.27 


Modal 

Frequencies 


1 


2 


3 


Experiment 

(Hz) 

Classical 
Theory  (Hz) 

Difference 

17.2 

15.8 

8.9% 

93.4 

94.4 

1.0% 

256.9 

280.9 

8.5% 

502.9 

524.4 

4.1% 

835.1 

901.7 

7.4% 

15.9 

15.6 

2.3% 

95.9 

90.4 

6.1% 

260.7 

273.1 

4.5% 

515.5 

508.4 

1.4% 

818.9 

828.6 

1.2% 

15.9 

15.3 

4.3% 

98.4 

86.8 

13.4% 

254.5 

254.4 

0.0% 

527.9 

494.4 

6.8% 

750.3 

724.3 

3.6% 

13.5 

15.0 

10.2% 

94.7 

83.8 

12.9% 

233.4 

228.9 

1.9% 

519.2 

486.5 

6.7% 

684.1 

661.1 

3.5% 

14.7 

14.4 

2.2% 

94.7 

81.3 

16.5% 

210.8 

203.1 

3.8% 

490.5 

465.4 

5.4% 

621.6 

638.8 

2.7% 

Figure  10.  Frequency  response  of  (a)  an  EFPI  fiber  optic  sensor  and  (b)  a  PZT  ceramic 
sensor  mounted  on  glass/epoxy  beam  containing  a  5.08  cm  delamination. 


Figure  14.  Frequency  response  of  an  EFPI  fiber  optic  sensor  mounted  on  a  glass/epoxy 
beam  containing  a  6.35  cm  delamination. 


C.  NEURAL  NETWORK  RESULTS 


The  neural  network  used  in  this  experiment  was  discussed  in  Section  IV.  Figure 
15  shows  the  variation  of  the  actual  modal  frequencies  used  to  test  the  network  for  the 
five  different  delamination  sizes:  1.27  cm,  2.54  cm,  3.81  cm,  5.08  cm,  and  6.35  cm.  Each 
delamination  was  incorporated  into  the  midplane  in  the  center  of  the  beam.  The  modal 
frequencies  are  normalized  to  a  maximum  of  100%.  Figure  16  shows  the  variation  of 
modal  frequencies  predicted  by  the  neural  network  for  the  five  different  delamination 
sizes.  Table  V  confirms  that  the  data  in  Figure  15  and  Figure  16  are  in  good  agreement. 
Table  V  compares  the  experimental  delamination  sizes  to  the  predicted  neural  network 
delamination  sizes.  The  average  error  of  the  neural  network  delamination  size  predictions 
was  5.9%.  The  error  values  ranged  from  0.9%  to  10.7%. 


Table  V:  Experimental  delamination  size  compared  to  predicted  neural  network 
delamination  size. 


Experimental 
Delamination  Size  (cm) 

Neural  Network 
Prediction  (cm) 

Percent 

Error 

1.27 

1.29 

1.2% 

2.54 

2.81 

10.6% 

3.81 

4.22 

10.7% 

5.08 

5.40 

6.3% 

6.35 

6.41 

0.9% 

The  neural  network  was  also  used  to  predict  delamination  location.  In  neural 
network  simulations,  four  randomly  selected  delamination  locations  for  each  of  the  five 
delamination  sizes  were  used  tested.  The  values  for  the  predicted  delamination  locations 
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are  shown  in  Table  VI.  The  average  error  between  the  actual  delamination  location  and 
the  neural  network  predicted  delamination  location  was  4.7%.  The  error  values  were 
evenly  distributed  ranging  from  zero  percent  to  13.9%. 


Table  VI:  Experimental  delamination  location  compared  to  predicted  neural  network 
delamination  location. 


Actual 


Actual 

Location  (cm) 

Predicted 
Location  (cm) 

Percent 

Error 

3.18 

3.20 

0.6% 

11.11 

10.92 

1.7% 

13.02 

11.76 

9.7% 

17.40 

17.65 

1.4% 

5.46 

6.22 

13.9% 

13.02 

13.32 

2.3% 

15.94 

16.55 

6.3% 

20.13 

19.42 

3.5% 

8.19 

9.12 

11.3% 

13.02 

13.74 

5.5% 

16.57 

17.15 

3.5% 

18.67 

19.33 

3.5% 

8.83 

8.83 

0.0% 

13.02 

13.54 

4.1% 

14.85 

15.55 

4.7% 

15.11 

15.07 

0.3% 

5.27 

5.74 

8.9% 

13.02 

12.52 

3.8% 

17.84 

17.88 

0.0% 

19.94 

19.71 

1.2% 
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Figure  17  shows  that  the  normalized  modal  frequencies  are  relatively  constant 
with  different  delamination  locations  for  a  2.54  cm  delamination  size.  Again,  Figure  18 
shows  that  the  neural  network  predicted  delamination  locations  are  in  good  agreement 
with  the  actual  delamination  location  for  a  2.54  cm  delamination  size  (see  figure  17). 
However,  figure  19  shows  that  the  normalized  modal  frequencies  change  with  different 
delamination  locations  for  a  6.35  cm  delamination.  Comparing  figure  19  and  figure  20 
confirms  that  the  neural  network  predicted  delamination  locations  agree  with  the  actual 
delamination  location  for  a  6.35  cm  delamination  size. 


Normalized  Frequency  (%) 


Normalized  Frequency  (%) 


Figure  16.  Neural  network  predicted  modal  frequencies  of  a  glass/epoxy  composite  beam 
for  five  prescribed  delamination  sizes. 


Normalized  Frequency  (%) 


5  7  9  11  13  15  17  19  21 

Actual  Delamination  Location  (cm) 


Figure  17.  Actual  modal  frequencies  of  four  different  2.54  cm  delamination  locations  in  a 
glass/epoxy  composite  beam. 


Normalized  Frequency  (%) 


Normalized 


Figure  19.  Actual  modal  frequencies  of  four  different  6.35  cm  delamination  locations  in  a 
glass/epoxy  composite  beam. 


Normalized  F 


Figure  20.  Neural  network  predicted  modal  frequencies  of  four  different  6.35  cm 
delamination  locations  in  a  glass/epoxy  composite  beam. 
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VI.  DISCUSSION 

A.  MODAL  FREQUENCIES  FOR  UNDAMAGED  COMPOSITE  BEAMS 

The  EFPI  fiber  optic  sensor  and  the  PZT  ceramic  sensor  gave  identical  modal 
frequencies  in  each  experiment  as  shown  in  Figures  8,  9,  and  10.  The  signal-to-noise 
ratio  was  higher  for  the  EFPI  fiber  optic  sensor  than  the  PZT  ceramic  patch,  however,  as 
the  butt  couple  changed  for  each  measurement,  the  signal-to-noise  ratio  varied  for  the 
EFPI  fiber  optic  sensor.  The  EFPI  fiber  optic  sensor  has  better  sensitivity  than  a  PZT 
ceramic  sensor  assuming  superior  fiber  optic  splices.  There  were  fine  structure 
differences  between  the  EFPI  fiber  optic  sensors  and  the  PZT  ceramic  sensors,  but  the 
major  structures  were  the  same. 

The  experimental  modal  frequencies,  calculated  classical  beam  theory  modal 
frequencies,  and  calculated  finite  element  model  frequencies  are  compared  in  Table  I. 
The  largest  percentage  difference  between  the  modal  frequencies  is  observed  in  the  first 
mode.  The  percent  difference  decreases  with  the  higher  order  modes.  The  percent 
difference  between  the  classical  beam  theory  and  the  finite  element  model  also  follows 
this  trend.  In  another  modal  frequency  study  the  opposite  was  observed  [31].  The  lower 
order  modes  had  the  least  percent  difference.  The  reason  for  this  difference  between 
studies  is  unknown. 

Beam  width  is  the  first  variable  used  to  experimentally  test  shifts  in  modal 
frequencies  predicted  by  the  classical  beam  theory  and  the  finite  element  model  used  in 
this  study.  Using  Table  II,  Figure  8,  and  Figure  9,  the  2.54  cm  wide  beam  clearly  has 
lower  modal  frequencies  than  the  2.33  cm  beam.  This  shift  is  expected  because  a  wider 
beam  is  stiffer  than  a  narrower  beam  assuming  identical  material  properties.  Both  the 
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classical  beam  theory  and  the  finite  element  model  predicted  that  modal  frequencies 
decrease  as  beam  width  increases.  However,  the  classical  beam  theory  and  the  finite 
element  model  over-predict  the  change  of  modal  frequencies  due  to  different  beam  widths 
as  seen  in  Table  II. 

B.  MODAL  FREQUENCIES  FOR  DELAMINATED  COMPOSITE  BEAMS 

Table  HI  compares  the  first  five  modal  frequeneies  of  an  undamaged  glass/epoxy 
composite  beam  to  five  other  beams  containing  various  sized  delaminations  incorporated 
in  their  midplanes.  All  the  beams  have  identical  material  properties  and  dimensions, 
since  variation  in  material  and  dimensions  can  cause  significant  modal  frequency  shifts  as 
shown  in  Table  H.  These  delamination-induced  shifts  can  be  seen  in  Figures  10  through 
14.  For  larger  delamination  sizes  and  high  modes,  the  modal  frequencies  are  down 
shifted  to  lower  frequencies.  As  delamination  size  increases,  the  modal  frequencies 
decrease. 

Table  IV  compares  the  modal  frequencies  obtained  experimentally  for  each 
delamination  size  to  the  modal  frequencies  calculated  by  the  classical  beam  theory.  The 
average  difference  is  5.6%  with  no  definite  pattern.  Therefore,  no  obvious  descripancies 
are  apparent  with  either  the  experimental  results  or  the  classical  beam  theory  calculations. 

C.  NEURAL  NETWORK  INSIGHTS  AND  IMPROVEMENTS 

The  normalized  predicted  neural  network  delamination  size  trends  shown  in 
Figure  16  are  very  similar  to  the  actual  normalized  delamination  size  trends  in  Figure  15. 
The  same  holds  true  for  the  actual  delamination  locations  (cf.  Figures  17  and  19) 
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compared  to  the  predicted  neural  network  delamination  locations  (cf.  Figures  18  and  20). 
These  matching  trends,  combined  with  a  5.9%  average  error  for  predicted  delamination 
sizes  and  a  4.7%  average  error  for  predicted  delamination  locations,  confirm  the  ability  of 
the  neural  network  to  accurately  predict  delamination  size  and  location. 

Figure  17  shows  that  the  normalized  modal  frequencies  do  not  change  with 
different  delamination  locations  for  a  2.54  cm  delamination  size.  Again,  Figure  18  shows 
that  the  neural  network  predicted  delamination  locations  are  in  good  agreement  with  the 
actual  delamination  location  for  a  2.54  cm  delamination  size.  However,  Figure  19  shows 
that  for  a  6.35  cm  delamination,  the  modal  frequencies  vary  with  different  delamination 
locations.  Higher  modes  show  the  most  variation.  Delamination  locations  farther  away 
from  the  clamped  end  of  the  beam  tend  to  vary  more  than  delaminations  close  to  the 
clamp.  Figure  20  again  confirms  that  the  neural  network  predicted  delamination 
locations  are  in  good  agreement  with  the  actual  delamination  locations.  Also,  increasing 
the  delamination  size,  from  2.54  cm  to  6.35  cm,  results  in  more  variation  of  normalized 
frequency  as  a  function  of  delamination  location.  By  comparing  Figures  15  through  20, 
the  neural  network  delamination  size  and  location  predictions  follow  the  actual 
delamination  sizes  and  locations  closely.  So,  apart  from  the  slight  oscillatory  behavior, 
the  neural  network  performed  successfully  in  predicting  the  prescribed  delamination  sizes 
and  delamination  locations  in  glass/epoxy  composite  beams. 
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VII.  CONCLUSIONS 

In  this  study,  EFPI  fiber  optic  sensors  were  used  to  measure  modal  frequencies  in 
composite  beams.  The  EFPI  fiber  optic  sensor  results  were  compared  to  results  from 
PZT  ceramic  sensors.  These  experimental  results  were  validated  by  classical  beam  theory 
and  a  finite  element  model.  A  feedforward  backpropagating  neural  network  was  used  to 
predict  the  experimental  results  based  on  the  classical  beam  theory  results. 

EFPI  fiber  optic  sensors  are  superior  to  the  PZT  ceramic  sensors  because  of 
sensitivity,  size,  ruggedness,  ability  to  be  embedded  into  the  composite  structure,  and 
immunity  to  electromagnetic  interference.  However,  the  sensitivity  is  a  function  of  the 
connection  of  the  EFPI  fiber  optic  sensor  to  the  fiber  optic  support  system.  The  EFPI 
fiber  optic  sensors  and  PZT  ceramic  sensors  measured  identical  modal  frequencies  in 
glass/epoxy  composite  beams. 

From  the  change  in  modal  frequencies  measured  by  the  EFPI  fiber  optic  sensors 
delamination  of  the  composite  beam  can  be  determined.  For  smaller  delamination  sizes 
the  first  five  modal  frequencies  increased  from  the  modal  frequencies  of  an  undamaged 
beam.  For  larger  delamination  sizes  the  first  two  modal  frequencies  increased  while  the 
next  three  modal  frequencies  decreased  compared  to  an  undamaged  beam.  Beam 
dimensions  also  have  an  effect  on  modal  frequencies.  Modal  frequencies  shift  up  as 
beam  width  or  thickness  is  increased. 

Classical  beam  theory  and  a  finite  element  model  were  used  to  compare  with  and 
validate  the  experimental  modal  frequencies.  Both  of  these  techniques  were  used  to 
validate  the  modal  frequencies  measured  by  the  EFPI  fiber  optic  sensors  and  the  PZT 
ceramic  sensors  for  the  undamaged  beams.  The  finite  element  model  did  not  account  for 
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beam  delaminations,  so  only  the  classical  beam  theory  could  be  used  to  compare  with  the 
experimental  modal  frequencies  from  the  damaged  beams. 

Using  modal  frequencies  calculated  by  MATLAB  simulations  based  on  classical 
beam  theory,  a  backpropagation  neural  network  was  trained  to  predict  delamination  size 
and  location  in  a  composite  beam.  The  error  of  the  neural  network  predicted 
delamination  sizes  ranged  from  0.9%  to  10.7%  with  an  average  error  of  5.9%.  The 
predicted  delamination  location  error  also  ranged  from  zero  percent  to  13.9%  with  an 
average  error  of  4.7%.  The  network  always  overpredicted  delamination  size  and 
generally  overpredicted  delamination  location.  However,  the  neural  network  successfully 
predicted  delamination  sizes  and  locations  in  glass/epoxy  composite  beams. 

The  neural  network  approach  is  well  suited  to  applications  that  predict  size  and 
location  of  delaminations  in  composite  lamina.  Further  work  related  to  this  study 
includes  improvements  to  the  accuracy  of  the  predictions  and  alternatives  to  the  specific 
neural  network  and  training  approaches.  Improvement  of  the  neural  network 
performance  can  be  achieved  by  training  the  network  with  more  information.  In 
particular,  incorporation  the  modal  shapes  are  possible.  The  modal  shapes  can  be  found 
experimentally,  with  the  classical  beam  theory,  and  with  the  finite  element  model.  Also 
improving  the  modeling  of  the  classical  beam  theory  equations  could  reduce  the  error  of 
the  calculated  modal  frequencies.  This  study  was  only  concerned  with  delaminations  in 
the  midplane  of  the  beam,  however,  this  technique  could  be  expanded  to  include 
prescribed  delaminations  located  between  other  laminate  layers.  The  neural  network 
architecture  could  also  be  changed,  e.g.,  number  of  layers  or  number  of  neurons. 


APPENDIX  A 


MATLAB  CODE  FOR  THE  CLASSICAL  BEAM  THEORY 
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The  following  appendix  contains  six  MATLAB  programs  and  one  text  file  that 
implement  the  classical  beam  theory  discussed  in  section  three  of  the  text.  The  first  program, 
Delam.m,  solves  16  coupled  equations.  The  only  input  into  the  program  is  beam  thickness.  The 
equations  are  solved  symbolically  so  they  can  be  used  for  any  application.  The  results  of 
Delam.m  can  be  seen  in  Function.txt.  The  second  program,  Property .m,  solves  for  the  material 
properties  of  the  beam.  Beam  width,  thickness,  length,  number  of  layers,  orientation  of  layers, 
and  the  lamina  properties  are  all  input  into  the  program  (all  parameters  are  in  English  units).  The 
material  properties  (Dl,  D2,  D3,  D4,  A2,  A3,  ml,  m2,  m3,  and  m4)  are  saved  into  a  workspace 
named  matrix.  Next,  Deter.m,  is  a  function  that  solves  for  the  determinate  of  the  16x16  matrix 
created  by  the  16  coupled  equations.  Deter.m  loads  the  workspace,  matrix,  and  uses  the  16 
equations  found  in  Function.txt  to  solve  the  determinate.  The  main  program  in  the  appendix  is 
Final.m.  This  program  uses  the  function  Deter.m  to  solve  for  five  roots,  i.e.,  modal  frequencies) 
of  1000  different  delaminations  at  up  to  11  different  locations  in  the  beam.  The  program  is  set¬ 
up  so  no  delamination  can  start  within  an  inch  of  the  clamp  holding  the  beam  or  end  more  than 
an  inch  from  the  free  end  of  the  beam.  The  program  takes  approximately  24  hours  to  run  on  a 
Pentium  133  MHz  processor.  A  1097  x  8  matrix  is  formed  and  saved  into  a  workspace  named 
final.  The  matrix  contains  the  first  five  modal  frequencies  associated  with  each  delamination, 
the  delamination  size,  and  the  beginning  and  ending  location  of  each  delamination.  Single.m  is  a 
program  identical  to  Final.m,  but  it  allows  the  user  to  solve  for  modal  frequencies  of  a  single 
delamination  size.  The  output  is  save  into  a  workspace  named  final2.  The  last  program  is 
Form.m.  This  program  changes  the  ascii  code  generated  by  MATLAB  into  a  text  file  that  can  be 
read  by  the  neural  network  software. 
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%Delam.m  Program 

%Modal  Frequency  Generation  for  Neural  Network 
%Generates  16  coupled  equations  symbolically 


%B21=B31; 

%B22=B32; 

%B23=B33; 

%B24=B34; 


%Segments  2  and  3  forced  to  vibrate  together 


%Beam  dimensions 
h=.061; 
h2=h/2; 
h3=h2; 


%lnitializing  constants 
B(1)=0; 

B(2)=0; 

B(3)=0; 

B(4)=0; 

B(5)=0; 

B(6)=0; 

B(7)=0; 

B(8)=0; 

B(9)=0; 

B(10)=0; 

B(11)=0; 

B(12)=0; 

B(13)=0; 

B(14)=0; 

B(15)=0; 

B(16)=0; 


%B1I 

%B12 

%B13 

%B14 

%B21,B31 

%B22,B32 

%B23,B33 

%B24,B34 

%B41 

%B42 

%B43 

%B44 

%C20 

%C21 

%C30 

%C31 


kl  =  sym('((ml*omega''2)/Dl)''(l/4)');  %Vibration  modes 
k2  =  sym('((m2*omega^2)/D2)^(l/4)'); 
k3  =  sym('((m3*omega''2)/D3)''(l/4)'); 
k4  =  sym('((m4*omega''2)/D4)''(l/4)'); 

u2  =  sym('B(13)+B(14)*x');  %Axial  modes 

u3  =  sym('B(15)+B(16)*x'); 


%Initial  equation 

W  =  sym('[B(l)*cosh(kl*x)+B(2)*sinh(kl*x)+B(3)*cos(kl*x)+B(4)*sin(kl*x); 
B(5)*cosh(k2*x)+B(6)*sinh(k2*x)+B(7)*cos(k2*x)+B(8)*sin(k2*x); 
B(5)*cosh(k3*x)+B(6)*sinh(k3*x)+B(7)*cos(k3*x)+B(8)*sin(k3*x); 
B(9)*cosh(k4*x)+B(10)*sinh(k4*x)+B(ll)*cos(k4*x)+B(12)*sin(k4*x)]'); 


%16  coupled  equations 

hi  =sym(W,l,l); 

h2  =  diff(hl,V); 

w4  =  sym(W,4,l); 

h3  =  diff(w4;x',2); 

h4  =  diff(w4,'x',3); 

h5  =  symsub(hl,sym(W,2,l)); 

dw2  =  diff(sym(W,2,l),'x'); 

h6  =  symsub(h2,dw2); 

VI  =  symmul('-Dr,diff(hl,'x',3)); 

V2  =  symmul('-D2',diff(sym(W,2,l),'x',3)); 
V3  =  symmul('-D3',diff(sym(W,3,l)/x',3)); 
V4  =  symmul('-D4',diff(w4,'x',3)); 

h7  =  symsub(Vl,symadd(V2,V3)); 


%Equation  1  @x=0 
%Equation  2 

%Equation  3  @x=L 
%Equation  4 
%Equation  5  @x=xl 

%Equation  6 

%Force  equations 


%Equation  7 


%Moment  Equations 


Ml  =  symmul('Dr,diff(hl,'x',2)); 

M2  =  symmul('D2',diff(sym(W,2,l),’x’,2)); 

M3  =  symmul('D3',diff(sym(W,3.1),'x',2)); 

M4  =  symmul('D4',diff(w4,'x',2)); 

P2  =  symmul('A2',diff(u2,'x’));  %Axial  Force  Equations 

P3  =  symmul('A3',diff(u3,'x')); 


h8a  =  symsub(Ml,symadd(M2,M3)); 

h8  =  synisub(h8a,symsub(symmul('(h/2-h3/2)',P3),symmuI('(h/2-h2/2)',P2))); 

%Equation  8 


h9  =  symadd(u2,symmul('(h/2-h2/2)',h2)); 
hlO  =  symsub(u3,syimnul('(h/2-h3/2)',h2)); 
hll  =  symadd(P2.P3); 


%Equation  9 

%Equation  10 
%Equation  11 


hl2  =  symsub(sym(W,2,l),w4);  ^Equation  12  @x=x2 

hl3  =  symsub(dw2,diff(w4,'x'));  %Equation  13 

hl4  =  symsub(V4.symadd(V2,V3));  %Equation  14 

hi  5a  =  symsub(M4,symadd(M2,M3)); 

hl5  =  symsub(hl5a,symsub(symmul('(h/2-h3/2)',P3),symmul('(h/2-h2/2)',P2))); 

%Equation  15 

hl6  =  symsub(u2,symsub(u3,symmul('(h/2)',diff(w4,'x')))); 

%Equation  16 


% - 

Function.txt 

Text  file  that  shows  the  16  coupled,  symbolic  equations  generated  by  Delam.m 
Equation  1 

hi  =B(l)*cosh(kl*x)+B(2)*sinh(kl*x)+B(3)*cos(kl*x)+B(4)*sin(kl*x) 

Equadon  2 

h2  =  B(l)*sinh(kl*x)*kl+B(2)»cosh(kl*x)*kl-B(3)*sin(kl*x)*kl+B(4)*cos(kl*x)*kl 


Equadon  3 

h3  =  B(9)*cosh(k4*x)*k4''2+B(10)*sinh(k4*x)*k4''2-B(l  l)*cos(k4*x)*k4''2- 
B(12)*sin(k4*x)*k4''2 


Equation  4 

h4  =  B(9)*sinh(k4*x)*k4''3+B(10)*cosh(k4*x)*k4''3+B(l  I)*sin(k4*x)*k4''3- 
B(12)*cos(k4*x)*k4''3 


Equation  5 

h5  =  B(l)*cosh(kl*x)+B(2)*sinh(kl*x)+B(3)*cos(kl*x)+B(4)*sin(kl*x)- 
B(5)*cosh(k2*x)-B(6)*sinh(k2*x)-B(7)*cos(k2*x)-B(8)*sin(k2*x) 

Equation  6 

h6  =  B(l)*sinh(kl*x)*kl+B(2)*cosh(kl*x)*kl-B(3)*sin(kl*x)*kl+B(4)*cos(kl*x)*kl- 
B(5)*sinh(k2*x)*k2-B(6)*cosh(k2*x)*k2+B(7)*sin(k2*x)*k2-B(8)*cos(k2*x)*k2 


Equation  7 

h7  =  -Dl*(B(l)*sinh(kl*x)*kl''3+B(2)*cosh(kl*x)*kl''3+B(3)*sin(kl*x)*kl''3- 
B(4)*cos(kl*x)*kl''3)+D2*(B(5)*sinh(k2*x)*k2''3+B(6)*cosh(k2*x)*k2''3 
+B(7)*sin(k2*x)*k2''3-B(8)*cos(k2*x)*k2''3)+D3*(B(5)*sinh(k3*x)*k3''3 
+B(6)*cosh(k3*x)*k3''3+B(7)*sin(k3*x)*k3''3-B(8)*cos(k3*x)*k3''3) 


Equation  8 

h8  =  Dl*(B(l)*cosh(kl*x)*kl''2+B(2)*sinh(kl*x)*kl''2-B(3)*cos(kl*x)*kl''2- 


B(4)*sin(kl*x)*kl''2)-D2*(B(5)*cosh(k2*x)*k2''2+B(6)*sinh(k2*x)*k2''2- 

B(7)*cos(k2*x)*k2''2-B(8)*sin(k2*x)*k2''2)-D3*(B(5)*cosh(k3*x)*k3''2 

+B(6)*sinh(k3*x)*k3''2-B(7)*cos(k3*x)*k3''2-B(8)*sin(k3*x)*k3''2)-(l/2*h- 

l/2*h3)*A3*B(16)+(l/2*h-l/2*h2)*A2*B(14) 


Equation  9 

h9  =  B(13)+B(14)*x+(l/2*h-l/2*h2)*(B(l)*sinh(kl*x)*kl+B(2)*cosh(kl*x)*kl- 
B(3)*sin(kl*x)*kl+B(4)*cos(kl*x)*kl) 


Equation  10 

hl0  =  B(15)+B(16)*x-(l/2*h-l/2*h3)*(B(l)*sinh(kl*x)*kl+B(2)*cosh(kl*x)*kl- 
B(3)*sin(kl*x)*kl+B(4)*cos(kl*x)*kl) 


Equation  11 

hll  =  A2*B(14)+A3*B(16) 


Equation  12 

hi  2  =  B(5)*cosh(k2*x)+B(6)*sinh(k2*x)+B(7)*cos(k2*x)+B(8)*sin(k2*x)- 
B(9)*cosh(k4*x)-B(10)*sinh(k4*x)-B(ll)*cos(k4*x)-B(12)*sin(k4*x) 


Equation  13 

hl3  =  B(5)*sinh(k2*x)*k2+B(6)*cosh(k2*x)*k2-B(7)*sin(k2*x)*k2 
+B(8)*cos(k2*x)*k2-B(9)*sinh(k4*x)*k4-B(10)*cosh(k4*x)*k4 
+B(ll)*sin(k4*x)*k4-B(12)*cos(k4*x)*k4 


Equation  14 

hl4  =  -D4*(B(9)*sinh(k4*x)*k4''3+B(10)’''cosh(k4*x)*k4''3+B(l  I)*sin(k4*x)*k4''3- 
B(12)*cos(k4»x)*k4''3)+D2*(B(5)*sinh(k2*x)*k2''3+B(6)*cosh(k2*x)*k2''3 
+B(7)*sin(k2*x)*k2''3-B(8)*cos(k2*x)*k2''3)+D3*(B(5)*sinh(k3*x)*k3''3 
+B(6)*cosh(k3*x)*k3''3+B(7)*sin(k3*x)*k3''3-B(8)*cos(k3*x)*k3^3) 


Equation  15 

hl5  =  D4*(B(9)*cosh(k4*x)*k4''2+B(10)*sinh(k4*x)*k4''2-B(l  l)*cos(k4*x)*k4''2- 
B(12)*sin(k4*x)*k4''2)-D2*(B(5)*cosh(k2*x)*k2''2+B(6)*sinh(k2*x)*k2''2- 
B(7)*cos(k2*x)*k2''2-B(8)*sin(k2*x)*k2''2)-D3*(B(5)*cosh(k3*x)*k3''2 
+B(6)*sinh(k3*x)*k3''2-B(7)*cos(k3*x)*k3''2-B(8)*sin(k3*x)*k3''2)-(l/2*h- 
l/2*h3)*A3*B(16)+(l/2*h-l/2*h2)*A2*B(14) 


Equation  16 

hl6  =  B(13)+B(14)*x-B(15)-B(16)*x+l/2*h*(B(9)*sinh(k4*x)*k4 

+B(10)’‘cosh(k4*x)*k4-B(ll)*sin(k4*x)*k4+B(12)*cos(k4*x)*k4) 


%  - . 

%Property.ra  Program 

%ModaI  Frequency  Generation  for  Neural  Network 
%Program  to  solve  for  the  material  properties  of  the  beam 


%lnputs 

b  = 

.917; 

h  = 

.061; 

L  = 

10.25; 

n  = 

8; 

%width 

%thickness 

%length 

%number  of  layers 


E 1 1  =  6. 14E6;  %  Lamina  properties 

E22=  1.7E6; 

G12=  1.03E6; 
vl2  =  .27; 
r  =  .18E-3; 


h2  =  h/2; 
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h3  =  h2; 
k  =  h/n; 
kk=h/(2*n); 


v21  =  (E22*vl2)/El  1 ;  %Constant  equations 

Qll  =Ell/(l-(vl2*v21)); 

Q12  =  (vl2*E22)/(l-(vl2*v21)); 

Q22  =  E22/(l-(vl2*v21)); 

Q66  =  G12; 

theta(l )  =  0;  %Layer  orientation 

theta(2)  =  pi/2; 

theta(3)  =  0; 

theta(4)  =  pi/2; 

theta(5)  =  pi/2; 

theta(6)  =  0; 

theta(7)  =  pi/2; 

theta(8)  =  0; 

%Loop  to  solve  Q  for  the  different  layer  orientations 
for  n  =  1:8 

Q(n)  =  Q1  l*(cos(theta(n)))''4+Q22*(sin(theta(n)))M  +2*(Q12+2*Q66)*(sin(theta(n)))''2*(cos(theta(n)))''2; 
end; 

%Solving  for  the  D's 

D=0; 

forj  =  1:8 

C  =  Q(j)*((((h/2)+(l-j)*k)''3)-((h/2)-j*k)^3); 

D  =  C+D; 
end; 

D1  =(1/3)*D; 

D4  =  D1; 

D=0; 
forj  =  1:4 

C  =  Q(j)*((((h/4)+(l-j)*kk)''3)-((h/4)-j*kk)''3); 

D  =  C+D; 
end; 

D2  =  (1/3)*D; 

D3  =  D2; 


%Solving  for  the  A's 
A  =  0; 
forj  =  1:4; 

C  =  Q(j)*(((h/4)+(l-j)*kk)-((h/4)-j*kk)); 

A  =  A+C; 

end; 


A2  =  1*A; 
A3  =  A2; 


%Solving  for  the  m's 
ml  =r*b*h; 
m4=ml; 
m2  =  r*b*h/2; 
m3=m2; 


save  matrix  D1  D2  D3  D4  A2  A3  ml  m2  m3  m4  L  h  h2  h3 
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% . . . . 

%Deter.m  Program 

%Modal  Frequency  Generation  for  Neural  Network 

%Function  to  determine  the  numeric  value  of  the  determinant  of  a  16x16  symbolic  matrix 

function  [d]  =  deter(y) 

load  matrix 
load  cor 


j=l:16; 

B=0*j; 


t=tm;  %Change  of  variables 

tt=ttm; 

omega  =  2*pi*y; 

kl  =  ((ml*omega''2)/Dl)''(l/4);  %Vibration  Modes 
k2  =  ((m2*omega''2)/D2)''(l/4); 
k3  =  ((m3*omega''2)/D3)''(l/4); 
k4  =  ((m4*omega‘'2)/D4)''(l/4); 

a(l,l)=l; 
a(1.3)=l; 
a(2,2)=l; 

a(2,4)=l; 
a(ll,14)=A2; 
a(ll,16)=A3; 


for  1=9:12 

B(i)=l; 

a(3,i)  =  B(9)*cosh(k4*x)*k4''2+B(10)*sinh(k4*x)*k4''2-B(ll)*cos(k4*x)*k4''2-B(12)*sin(k4*x)*k4''2: 
a(4,i)  =  B(9)*sinh(k4*x)*k4''3+B(10)*cosh(k4*x)*k4''3+B(l  l)*sin(k4*x)*k4^3-B(12)*cos(k4*x)*k4''3; 
B(i)=0; 


%Initializes  matrix  values 


%Starts  loop  for  x=L 


x=tt; 

for  1=9: 1 6  %Starts  loop  for  x  equal  to  end  of  delamination 

B(i)=l; 

a(16,i)  =  B(13)+B(14)*x-B(15)-B(16)*x+(l/2*h)*(B(9)*sinh(k4*x)*k4+B(10)*cosh(k4*x)*k4- 
B(ll)*sin(k4*x)*k4+B(12)*cos(k4*x)*k4); 

B(i)=0; 

end; 

x=t; 

for  1=1 :8  %Starts  loop  for  x  equal  to  front  of  delamination 

B(i)=l; 

a(5,i)  =  B(l)*cosh(kl*x)+B(2)*sinh(kl*x)+B(3)*cos(kl*x)+B(4)*sin(kl*x)-B(5)*cosh(k2*x)-B(6)*smh(k2*x)- 

B(7)*cos(k2*x)-B(8)*sin(k2*x); 

a(6,i)  =  B(l)*sinh(kl*x)*kl+B(2)*cosh(kl*x)*kl-B(3)*sin(kl*x)*kl+B(4)*cos(kl*x)*kl-B(5)*sinh(k2*x)*k2- 
B(6)*cosh(k2*x)*k2+B(7)*sin(k2*x)*k2-B(8)*cos(k2*x)*k2; 

a(7,i)  =  -Dl*(B(l)*sinh(kl*x)*kl''3+B(2)*cosh(kl*x)*kl''3+B(3)*sin(kl*x)*kl''3- 
B(4)*cos(kl*x)*kl''3)+D2*(B(5)*sinh(k2*x)*k2''3+B(6)*cosh(k2*x)*k2''3+B(7)*sin(k2*x)*k2''3- 

B(8)*cos(k2*x)*k2''3)+D3*(B(5)*sinh(k3*x)*k3''3+B(6)*cosh(k3*x)*k3''3+B(7)*sin(k3*x)*k3''3-B(8)*cos(k3*x)*k3''3); 

B(i)=0; 

end; 

for  i=l:16 

B(i)=l; 

x=t; 

a(8,i)  =  Dl*(B(l)*cosh(kl*x)*kl''2+B(2)*sinh(kl*x)*kl''2-B(3)*cos{kl*x)*kl''2-B(4)*sin(kl*x)*kl''2)- 
D2*(B(5)*cosh(k2*x)*k2''2+B(6)*sinh(k2*x)*k2''2-B(7)*cos(k2*x)*k2''2-B(8)*sin(k2*x)*k2''2)- 
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D3*(B(5)*cosh(k3*x)*k3''2+B(6)*sinh(k3*x)*k3''2-B(7)*cos(k3*x)*k3''2-B(8)*sin(k3*x)*k3''2)-(l/2*h- 

l/2*h3)*A3*B(16)+(l/2*h-l/2*h2)*A2*B(I4); 

a(9.i)  =  B(13)+B(14)*x+(l/2*h-l/2*h2)*(B(l)*sinh(kl*x)*kl+B(2)*cosh{kl*x)*kl-B(3)*sin(kl*x)*kl+B(4)*cos(kl*x)*kl); 
a(10,i)  =  B(15)+B(16)*x-(l/2*h-l/2*h3)*(B(l)*sinh(kl*x)*kl+B(2)*cosh(kl*x)*kl-B(3)*sin(kl*x)*kl+B(4)*cos(kl*x)*kl); 
x=tt; 

a(15,i)  =  D4*(B(9)*cosh(k4*x)*k4''2+B(10)*sinh(k4*x)*k4'^2-B(ll)*cos(k4*x)*k4''2-B(12)*sin(k4*x)*k4''2)- 
D2*(B(5)*cosh(k2*x)*k2''2+B(6)*sinh(k2*x)*k2''2-B(7)*cos(k2*x)*k2''2-B(8)*sin(k2*x)*k2''2)- 
D3*(B(5)*cosh(k3*x)*k3''2+B(6)*sinh(k3*x)*k3''2-B(7)*cos(k3*x)*k3''2-B(8)*sin(k3*x)*k3''2)-(l/2*h- 
l/2*h3)*A3*B(16)+(l/2*h-l/2*h2)*A2*B(14); 

B(i)=0; 

end; 


x=tt; 

for  i=5:12 
B(i)=l; 

a(12,i)  =  B(5)*cosh(k2*x)+B(6)*sinh(k2*x)+B(7)*cos(k2*x)+B(8)*sin(k2*x)-B(9)*cosh(k4*x)-B(10)*sinh(k4*x)- 
B(ll)*cos(k4*x)-B(12)*sin(k4*x); 

a(13,i)  =  B(5)*sinh(k2*x)*k2+B(6)*cosh(k2*x)*k2-B(7)*sin(k2*x)*k2+B(8)*cos(k2*x)*k2-B(9)*sinh(k4*x)*k4- 

B(10)*cosh(k4*x)*k4+B(ll)*sin(k4*x)*k4-B(12)*cos(k4*x)*k4; 

a(14,i)  =  -D4*(B(9)*sinh(k4*x)*k4''3+B(10)*cosh(k4*x)*k4''3+B(l  I)*sin(k4*x)*k4''3- 

B(12)*cos(k4*x)*k4''3)+D2*(B(5)*sinh(k2*x)*k2''3+B(6)*cosh(k2*x)*k2''3+B(7)*sin(k2*x)*k2^3- 

B(8)*cos(k2*x)*k2''3)+D3*(B(5)*sinh(k3*x)*k3''3+B(6)*cosh(k3*x)*k3'‘3+B(7)*sin(k3*x)*k3''3-B(8)*cos(k3*x)*k3''3); 

B(i)=0; 

end; 

d  =det(a); 


% - 

%Final.m  Program 

%Modal  Frequency  Generation  for  Neural  Network 


%Finds  the  first  five  modal  frequencies  for  different  delamination  sizes  and  locations  in  a  composite  beam 
freq  =  [0,0,0,0,0,0,0,0,0,0,0,0,0.0,0,0,0,0,0,0,0,0,01;  %  watch  out  for  the  size  of  these 

matrices 

uuuu=  [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01 ; 
uuu=  [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01 ; 
uu  =  [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01; 

u  =  [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01; 


range  =  1000; 
nmaxl  =  20; 
nmax2  =10; 
nmax3  =  10; 


%  the  frequencies  are  between  0  and  range 
%  number  of  initial  intervals 

%  range/(nmaxl*nmax2*nmax3*nmax4)  is  the  accuracy  on  the  frequency 
%  make  sure  you  put  round  numbers  otherwise  you  get  "truncated  errors" 
from  maple 


nmax4  =10; 
p  =  nmaxl/range; 

moots  =  5;  %  number  of  roots  to  find 


%Loop  to  calculate  1000  different  delaminations  and  locations 
Z=ll; 

C=0; 
q=0; 

L=10,25; 
x(2)=l; 

x(3)=8.25/(Z-l); 
x(4)=2*(8.25/(Z-l)); 
x(5)=3*(8.25/(Z-l)); 
x(6)=4*(8.25/(Z-l)); 
x(7)=5*(8.25/(Z-l)); 
x(8)=6*(8.25/(Z-l)); 
x(9)=7*(8,25/(Z-l)); 


%Delamination  locations 
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x(10)=8*(8.25/(Z-l)); 

x(ll)=9*(8.25/(Z-l)); 

for  g=l  :200  %Changes  delamination  size  by  0.04'' 

t  =  5.145-(.02*g); 
tt  =  5.10501+(.02'^g); 

delam=tt-t;  %Delamination  size 

x(l)=t; 

if  q<Z 

C=C+1; 

end; 


for  q=l  :Z;  %Changes  delamination  location 

tm=x(q); 

ttm=delam+x(q); 

save  cor  tm  ttm 

if  ttm  <  9.25  %If  delamination  extends  out  of  beam 


% . 

%Prgram  to  find  n  (numeric)  roots  of  the  determinant  of  16x16  matrix  (symbolic) 
i=l; 

m=l; 

u(l)  =  deter(O.OOl); 
for  j=l:nmaxl 

u(j+l)  =  deter(j/p); 

end; 

while  m  <  nmaxl+1 


index=0; 

if  u(m)'*u(m+l)  <  0 

uu(l)  =  u(m); 
uu(nmax2)  =  u(m+l); 
k=l; 

while  k<  nmax2 

y  =  (m-l+k/nmax2)/p; 

ifk<nmax2-l 

uu(k+l)=deter(y); 

end; 

if  uu(k)'''uu(k+l)  <  0 

uuu(l)=uu(k); 

uuu(nmax3)=uu(k+l); 

1=1; 

while  1  <  nmax3 

yy  =  (m-l+(k-l+ 

(l/nmax3))/nmax2)/p; 
if  1  <  nmax3-l 

uuu(l+ 1  )=deter(y  y) ; 

end; 


if  uuu(l)'''uuu(l+l)  <  0 

uuuu(l)  =  uuu(l); 
uuuu(nmax4)  =  uuu(l+l); 
11=1; 

while  11  <  nmax4 

yyy  =  (m-l+(k-l+ 

((ll)+ll/nmax4)/ 

nmax3)/nmax2)/p; 


deter(yyy); 


if  11  <  nmax4-l 
uuuu(ll+l)  = 

end; 

if  uuuu(ll)'*'uuuu(ll+l)  <  0; 

freq(i)  =  (2*yyy 


-  l/(nmax2'''nmax3 


‘nmax4*p))/2; 


i=i+l; 

index=l; 

break; 


end; 

if  index  ==  I 

break; 

end; 

1=1+1; 

end; 

end; 

if  index  ==  1 
break; 
end; 
k=k+l; 

end; 

if  i  ==  nroots+1 
break; 

end; 

end; 


m=m+l; 

end; 

% — . 

MF(C,l)=freq(l); 

%Loads  values  into  matrix  MF 

MF(C,2)=freq(2); 

MF(C,3)=freq(3); 

MF(C,4)=freq(4); 

MF(C,5)=freq(5); 

MF(C,6)=delam; 

MF(C,7)=tm; 

MF(C,8)=ttm; 

C=C+1; 

end; 

%Ends  if  loop 

end; 

%Ends  location  loop 

end; 

%Ends  size  loop 

S  =  size(MF); 

save  final  MF  S 

%Saves  workspace  into  final 

end; 

11  =  11+1; 

end; 


% . . . 

Single.ra  Program 

%Modal  Frequency  Generation  for  Neural  Network 

%Does  the  same  operation  as  final.m,  however  allows  for  the  determination  of  modal  frequencies  for  single  delamination 


freq  =  [0,0, 0,0, 0,0, 0,0, 0,0,0, 0,0, 0,0, 0,0, 0,0,0, 0,0,0]; 

uuuu=  [0,0,0.0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] ; 
uuu=  [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01 ; 
uu  =  [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01; 
u  =  [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01; 


%  watch  out  for  the  size  of  these 
matrices 


range  =  1000; 
nmaxl  =  20; 
nmax2  =10; 
nmax3  =  10; 
nmax4  =  10; 
p  =  nmaxl/range; 
moots  =  5; 


I  the  frequencies  are  between  0  and  range 
I  number  of  initial  intervals 

)  range/(nmaxl*nmax2*nmax3*nmax4)  is  the  accuracy  on  the  frequency 
)  make  sure  you  put  round  numbers  otherwise  you  get  "truncated  errors"  from  maple 


)  number  of  roots  to  find 
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%Loop  to  calculate  1000  different  delaminations  and  locations 
Z=ll; 

C=l; 

q=0; 

L=10.25; 

x(2)=l;  %Delamination  locations 

x(3)=8.25/(Z-l); 

x(4)=2*(8.25/(Z-l)); 

x(5)=3*(8.25/(Z-l)); 

x(6)=4*(8.25/(Z-l)); 

x(7)=5*(8.25/(Z-l)); 

x(8)=6*(8.25/(Z-l)); 

x(9)=7*(8.25/(Z-l)); 

x(10)=8*(8.25/(Z-l)); 

x(ll)=9*(8.25/(Z-l)); 

%Enter  delamination  size  below 

t  =  5.125-(0.5/2);  %Changes  delamination  size  by  0.1” 

tt  =  5.125+(0.5/2); 

delam=tt-t;  %Delamination  size 

x(l)=t; 


for  q=l  :Z;  %Changes  delamination  location 

%if  q<Z 

%C=C+1; 

%end; 

tm=x(q); 
ttm=delam+x(q); 
save  cor  tm  ttm 
if  ttm  >  9.25 
C=C+1; 
end; 

if  ttm  <  9.25  %If  delamination  extends  out  of  beam 


%prgram  to  find  n  (numeric)  roots  of  the  determinant  of  16x16  matrix  (symbolic) 

i=l; 

m=l; 

u(l)  =  deter(O.OOl); 
for  j=l:nmaxl 

u(j+l)  =  deter(j/p); 

end; 

while  m  <  nmaxl+1 
index=0; 

if  u(m)*u(m+l)  <  0 

uu(l)  =  u(m); 
uu(nmax2)  =  u(m+l); 
k=l; 

while  k<  nmax2 

y  =  (m-l+k/nmax2)/p; 

if  k  <  nmax2-l 

uu(k+l)=deter(y); 

end; 

if  uu(k)*uu(k+l)  <  0 

uuu(l)=uu(k); 

uuu(nmax3)=uu(k+l); 

1=1; 

while  1  <  nmax3 

yy  =  (m-l+(k-l+ 


(l/nmax3))/nmax2)/p; 


ifi== 

end; 

end; 

m=m+l; 

end; 

% - 

MF(C,l)=freq(l) 

MF(C,2)=freq(2) 

MF(C,3)=freq(3) 

MF(C,4)=freq(4) 

MF(C,5)=freq(5) 

MF(C,6)=delam; 

MF(C,7)=tm; 

MF(C,8)=ttm; 

C=C+1; 

end; 

end; 

%end; 

S  =  size(MF); 
save  final  MF  S 


if  I  <  nmax3-l 

uuu(l+l  )=deter(yy ) ; 

end; 

if  uuu(l)*uuu(l+l)  <  0 

uuuu(l)  =  uuu(l); 
uuuu(nmax4)  =  uuu(l+l); 

11=1; 

while  11  <  nmax4 

yyy  =  (m-l+(k-l+ 

((l-l)+ll/nmax4)/ 

nmax3)/nmax2)/p; 


end; 

nroots+l 

break; 


end; 

if  index  ==  1 

break; 

end; 

k=k+l; 


if  ll<nmax4-l 
uuuu(ll+l)  = 

deter(yyy); 

end; 

if  uuuu(ll)*uuuu(ll+l)  <  0; 

freq(i)  =  (2*yyy 

-l/(nmax2*nmax3 

*nmax4*p))/2; 

i=i+l; 

index=l; 

break; 

end; 

11  =  11+1; 
end; 

end; 

if  index  ==  1 


break; 

end; 

1=1+1; 

end; 


%loads  values  into  matrix  MF 


%ends  if  loop 
%ends  location  loop 
%ends  size  loop 

%saves  workspace  into  final 
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% - - 

%Form.m  Program 

%Modal  Frequency  Generation  for  Neural  Network 

%Changes  ascii  code  to  text  file  so  the  five  modal  frequencies,  delamination  size,  and  delamination  location  can  be  used  to  train  the 
neural  network  (from  average.mat) 

fod  =  fopen('E:\FARHAD\GIL\IMPORTAN\outtest2.txt',’w+'); 
load  average 

outcount  =  fprintf(fod,'%6.3f  %6.3f  %6.3f  %6.3f  %6.3f  %6.3f  %6.3f\n',MN') 
fclose(fod); 


APPENDIX  B 

MODAL  FREQUENCIES  GENERATED  BY  THE  CLASSICAL  BEAM 


THEORY 
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The  following  appendix  contains  the  data  generated  using  the  MATLAB  codes 
implementing  the  classical  beam  theory.  The  first  five  columns  are  the  first  five  modal 
frequency  for  a  glass-epoxy  composite  beam.  The  units  on  modal  frequencies  are  hertz 
(Hz).  The  sixth  column  is  the  delamination  size.  The  delamination  size  is  measured  in 
inches  (the  rest  of  the  thesis  uses  MKS  units).  The  final  column  is  the  delamination 
location,  also  measured  in  inches.  The  location  is  described  at  the  center  of  the 
delamination.  It  was  found  by  averaging  the  beginning  and  ending  of  each  delamination 


location. 
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Mode  ] 

16.075 

16.075 

16.075 

16.075 

16.075 

16.075 

16.075 

16.075 

16.075 

16.075 

16.075 

16.075 

15.925 

15.925 

15.925 

16.025 

16.025 

16.025 

16.075 

16.075 

16.075 

16.075 

16.025 

15.825 

15.825 

15.875 

15.925 

15.925 

16.025 

16.025 

16.075 

16.075 

16.075 

16.025 

15.725 

15.725 

15.775 

15.875 

15.925 

15.925 

16.025 

16.025 

16.075 

16.075 

16.025 

15.625 

15.575 

15.725 

15.825 

15.875 

15.925 

15.925 

16.025 

16.075 

16.075 

15.925 

15.525 

15.425 

15.625 

15.725 

15.825 


DATA  GENETERATED  USING  CLASSICAL  BEAM  THEORY 


Mode  2 

Mode  3 

Mode  4 

Mode  5 

Size 

Location 

100.775 

282.125 

552.875 

913.925 

0.000 

5.125 

100.775 

282.125 

552.875 

913.925 

0.000 

1.000 

100.775 

282.125 

552.875 

913.925 

0.000 

0.825 

100.775 

282.125 

552.875 

913.925 

0.000 

1.650 

100.775 

282.125 

552.875 

913.925 

0.000 

2.475 

100.775 

282.125 

552.875 

913.925 

0.000 

3.300 

100.775 

282.125 

552.875 

913.925 

0.000 

4.125 

100.775 

282.125 

552.875 

913.925 

0.000 

4.950 

100.775 

282.125 

552.875 

913.925 

0.000 

5.775 

100.775 

282.125 

552.875 

744.025 

0.000 

6.600 

100.775 

282.125 

650.075 

755.075 

0.000 

7.425 

100.225 

282.125 

550.075 

913.925 

0.040 

5.105 

100.425 

281.925 

552.825 

913.125 

0.040 

1.000 

100.375 

281.775 

552.775 

913.825 

0.040 

0.825 

100.725 

282.025 

551.375 

909.425 

0.040 

1.650 

100.725 

281.225 

550.625 

913.025 

0.040 

2.475 

100.625 

280.925 

552.525 

912.175 

0.040 

3.300 

100.375 

281.525 

552.175 

909.425 

0.040 

4.125 

100.275 

282.125 

550.175 

913.675 

0.040 

4.950 

100.225 

281.725 

551.325 

911.025 

0.040 

5.775 

100.375 

280.825 

552.825 

910.425 

0.040 

6.600 

100.525 

280.525 

550.925 

915.025 

0.040 

7.425 

94.425 

282.125 

544.425 

913.875 

0.080 

5.085 

100.225 

281.875 

552.775 

912.175 

0.080 

1.000 

100.075 

281.425 

552.725 

913.575 

0.080 

0.825 

100.675 

281.875 

544.425 

906.225 

0.080 

1.650 

100.725 

280.325 

544.425 

912.325 

0.080 

2.475 

100.425 

279.425 

552.275 

910.225 

0.080 

3.300 

100.025 

280.925 

551.325 

905.525 

0.080 

4.125 

94.425 

282.125 

544.425 

913.575 

0.080 

4.950 

94.425 

281.325 

544.425 

908.025 

0.080 

5.775 

94.425 

279.425 

552.825 

907.175 

0.080 

6.600 

100.325 

278.925 

544.425 

910.225 

0.080 

7.425 

94.425 

282.125 

544.425 

913.725 

0.120 

5.065 

100.025 

281.775 

552.675 

910.925 

0.120 

1.000 

94.425 

281.225 

552.675 

913.125 

0.120 

0.825 

100.625 

281.725 

544.425 

902.575 

0.120 

1.650 

100.725 

279.375 

544.425 

911.775 

0.120 

2.475 

100.275 

278.675 

552.075 

908.075 

0.120 

3.300 

94.425 

280.525 

550.425 

901.825 

0.120 

4.125 

94.425 

282.075 

544.425 

913.425 

0.120 

4.950 

94.425 

280.875 

544.425 

904.425 

0.120 

5.775 

94.425 

278.225 

552.725 

904.325 

0.120 

6.600 

100.125 

277.425 

544.425 

914.075 

0.120 

7.425 

94.425 

282.075 

542.375 

913.525 

0.160 

5.045 

94.425 

281.675 

552.425 

909.375 

0.160 

1.000 

94.425 

281.025 

552.575 

912.375 

0.160 

0.825 

100.575 

281.525 

544.425 

894.425 

0.160 

1.650 

100.675 

278.425 

544.425 

911.275 

0.160 

2.475 

100.075 

277.575 

551.925 

905.625 

0.160 

3.300 

94.425 

280.075 

544.425 

894.425 

0.160 

4.125 

94.425 

282.075 

542.425 

913.225 

0.160 

4.950 

94.425 

280.375 

544.425 

901.775 

0.160 

5.775 

94.425 

276.925 

552.575 

901.675 

0.160 

6.600 

94.425 

275.925 

544.425 

912.625 

0.160 

7.425 

94.425 

282.025 

540.025 

913.125 

0.200 

5.025 

94.425 

281.625 

552.125 

907.575 

0.200 

1.000 

94.425 

280.825 

552.375 

911.375 

0.200 

0.825 

100.575 

281.275 

544.425 

894.425 

0.200 

1.650 

100.625 

277.425 

543.175 

910.825 

0.200 

2.475 

94.425 

276.575 

551.725 

903.025 

0.200 

3.300 
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Mode  1 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Size 

Location 

15.875 

91.875 

265.275 

526.325 

852.625 

0.880 

5.775 

15.925 

94.225 

255.925 

524.425 

852.075 

0.880 

6.600 

16.025 

94.425 

256.375 

507.775 

840.425 

0.880 

7.425 

15.625 

91.075 

275.025 

510.425 

844.425 

0.920 

4.665 

14.075 

94.425 

272.225 

520.425 

858.725 

0.920 

1.000 

13.925 

94.425 

271.525 

520.375 

854.425 

0.920 

0.825 

14.425 

94.425 

268.575 

516.775 

855.825 

0.920 

1.650 

14.425 

94.425 

261.875 

522.625 

844.425 

0.920 

2.475 

15.175 

94.425 

264.425 

526.675 

844.425 

0.920 

3.300 

15.425 

92.175 

273.075 

514.425 

844.425 

0.920 

4.125 

15.725 

90.825 

273.875 

511.875 

844.425 

0.920 

4.950 

15.875 

91.575 

264.175 

523.675 

851.125 

0.920 

5.775 

15.925 

94.025 

254.425 

521.425 

844.425 

0.920 

6.600 

16.025 

94.425 

255.575 

505.925 

835.925 

0.920 

7.425 

15.625 

90.725 

274.125 

509.375 

836.925 

0.960 

4.645 

14.025 

94.425 

271.025 

517.925 

857.275 

0.960 

1.000 

13.925 

94.425 

270.325 

517.425 

852.425 

0.960 

0.825 

14.375 

94.425 

267.525 

515.825 

851.625 

0.960 

1.650 

14.425 

94.425 

261.225 

520.725 

839.175 

0.960 

2.475 

15.175 

94.425 

264.425 

523.325 

844.425 

0.960 

3.300 

15.425 

91.775 

272.425 

513.025 

841.275 

0.960 

4.125 

15.725 

90.425 

272.775 

510.825 

838.825 

0.960 

4.950 

15.875 

91.275 

263.075 

520.825 

844.425 

0.960 

5.775 

15.925 

93.775 

253.925 

518.175 

839.425 

0.960 

6.600 

16.025 

94.425 

254.425 

504.075 

830.775 

0.960 

7.425 

15.625 

90.375 

273.125 

508.375 

828.625 

1.000 

4.625 

13.925 

94.425 

269.425 

515.625 

856.075 

1.000 

1.000 

13.825 

94.425 

269.025 

514.425 

850.425 

1.000 

0.825 

14.325 

94.425 

266.425 

514.425 

844.425 

1.000 

1.650 

14.425 

94.425 

260.575 

518.625 

831.825 

1.000 

2.475 

15.125 

94.425 

263.925 

519.425 

844.425 

1.000 

3.300 

15.425 

91.425 

271.625 

511.125 

835.025 

1.000 

4.125 

15.675 

90.125 

271.625 

509.425 

831.875 

1.000 

4.950 

15.875 

90.925 

261.925 

517.675 

844.425 

1.000 

5.775 

15.925 

93.575 

253.025 

514.425 

832.425 

1.000 

6.600 

16.025 

94.425 

254.125 

502.275 

830.025 

1.000 

7.425 

15.575 

90.075 

272.025 

507.375 

820.175 

1.040 

4.605 

13.875 

94.425 

268.425 

513.375 

855.025 

1.040 

1.000 

13.775 

94.425 

267.675 

511.675 

844.425 

1.040 

0.825 

14.275 

94.425 

265.375 

514.175 

841.125 

1.040 

1.650 

14.425 

94.425 

259.425 

516.275 

824.425 

1.040 

2,475 

15.075 

94.125 

263.425 

516.175 

844.425 

1.040 

3.300 

15.425 

91.075 

270.825 

509.325 

828.225 

1.040 

4.125 

15.675 

89.425 

270.425 

508.525 

824.425 

1.040 

4,950 

15.875 

90.675 

260.775 

514.275 

844.425 

1.040 

5.775 

15.925 

93.375 

252.125 

511.175 

825.175 

1.040 

6.600 

16.025 

94.425 

253.375 

500.525 

820.025 

1.040 

7.425 

15.575 

89.425 

270.875 

506.425 

811.425 

1.080 

4,585 

13.825 

94.425 

267.025 

511.175 

854.075 

1.080 

1.000 

13.725 

94.425 

266.275 

508.925 

844.425 

1.080 

0.825 

14.225 

94.425 

264.325 

513.425 

834.425 

1.080 

1.650 

14.425 

94.425 

259.425 

513.675 

817.425 

1.080 

2.475 

15.075 

93.775 

262.925 

512.375 

844.175 

1.080 

3.300 

15.425 

90.725 

269.425 

507.625 

820.925 

1.080 

4.125 

15.675 

89.425 

269.175 

507.225 

818.125 

1.080 

4.950 

15.875 

90.425 

259.425 

510.675 

844.425 

1.080 

5.775 

15.925 

93.175 

251.225 

507.525 

817.525 

1.080 

6.600 

16.025 

94.425 

252.675 

494.425 

810.025 

1.080 

7.425 

15.525 

89.425 

269.425 

505.525 

802.725 

1.120 

4.565 

13.775 

94.425 

265.625 

509.175 

853.175 

1.120 

1.000 

13.675 

94.425 

264.425 

506.325 

844.425 

1.120 

0.825 

14.175 

94.425 

263.275 

512.675 

827.825 

1.120 

1.650 

14.425 

94.425 

258.875 

510.775 

810.625 

1.120 

2.475 

Mode  1 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Size 

Location 

15.025 

93.425 

262.425 

508.425 

843.025 

1.120 

3.300 

15.375 

90.375 

268.925 

506.025 

813.175 

1.120 

4.125 

15.675 

89.175 

267.825 

505.875 

811.425 

1.120 

4.950 

15.875 

90.175 

258.425 

506.775 

844.425 

1.120 

5.775 

15.925 

93.025 

250.375 

503.725 

809.425 

1.120 

6.600 

16.025 

94.425 

252.025 

494.425 

800.025 

1.120 

7.425 

15.525 

89.125 

268.375 

504.425 

793.875 

1.160 

4.545 

13.725 

94.425 

264.125 

507.275 

852.325 

1.160 

1.000 

13.625 

94.425 

263.275 

503.825 

844.425 

1.160 

0.825 

14.125 

94.425 

262.275 

511.925 

820.175 

1.160 

1.650 

14.425 

94.425 

258.375 

507.625 

804.175 

1.160 

2.475 

15.025 

93.075 

261.875 

504.425 

841.875 

1.160 

3.300 

15.375 

90.025 

267.825 

504.425 

805.075 

1.160 

4.125 

15.625 

88.825 

266.425 

504.425 

805.075 

1.160 

4.950 

15.825 

89.425 

257.275 

502.675 

843.075 

1.160 

5.775 

15.925 

92.825 

244.425 

494.425 

801.625 

1.160 

6.600 

16.025 

94.425 

251.325 

494.425 

792.625 

1.160 

7.425 

15.425 

88.825 

267.025 

503.825 

785.075 

1.200 

4.525 

13.675 

94.425 

262.625 

505.525 

851.425 

1.200 

1.000 

13.575 

94.425 

261.675 

501.425 

844.075 

1.200 

0.825 

14.075 

94.425 

261.225 

511.175 

811.875 

1.200 

1.650 

14.425 

94.425 

257.875 

504.175 

794.425 

1.200 

2.475 

14.425 

92.725 

261.325 

500.425 

840.725 

1.200 

3.300 

15.325 

89.425 

266.675 

503.075 

794.425 

1.200 

4.125 

15.625 

88.525 

264.425 

502.825 

794.425 

1.200 

4.950 

15.825 

89.425 

256.075 

494.425 

841.425 

1.200 

5.775 

15.925 

92.625 

244.425 

494.425 

793.425 

1.200 

6.600 

16.025 

94.425 

250.675 

494.175 

784.275 

1.200 

7.425 

15.425 

88.525 

265.625 

503.025 

776.325 

1.240 

4.505 

13.625 

94.425 

261.075 

503.925 

850.425 

1.240 

1.000 

13.525 

94.425 

260.075 

494.425 

843.325 

1.240 

0.825 

14.025 

94.425 

260.225 

510.375 

803.175 

1.240 

1.650 

14.425 

94.425 

257.425 

500.425 

792.575 

1.240 

2.475 

14.425 

92.375 

260.675 

494.425 

839.425 

1.240 

3.300 

15.325 

89.325 

265.425 

501.725 

788.075 

1.240 

4.125 

15.625 

88.225 

263.425 

501.075 

793.275 

1.240 

4.950 

15.825 

89.375 

254.425 

493.925 

839.425 

1.240 

5.775 

15.925 

92.425 

244.425 

491.875 

785.275 

1.240 

6.600 

16.025 

94.425 

250.025 

492.725 

777.275 

1.240 

7.425 

15.425 

88.275 

264.125 

502.225 

767.725 

1.280 

4.485 

13.575 

94.425 

259.425 

502.425 

844.425 

1.280 

1.000 

13.425 

94.425 

258.425 

494.425 

842.675 

1.280 

0.825 

13.925 

94.425 

259.225 

509.425 

794.025 

1.280 

1.650 

14.425 

94.425 

256.925 

494.425 

787.525 

1.280 

2.475 

14.425 

92.025 

260.025 

492.425 

838.125 

1.280 

3.300 

15.325 

89.025 

264.075 

500.525 

779.375 

1.280 

4,125 

15.625 

87.925 

261.875 

494.425 

787.925 

1.280 

4.950 

15.825 

89.125 

253.675 

489.275 

837.575 

1.280 

5.775 

15.925 

92.275 

244.425 

487.825 

777.075 

1.280 

6.600 

16.025 

94.425 

244.425 

491.325 

769.425 

1.280 

7.425 

15.425 

87.925 

262.525 

501.425 

759.275 

1.320 

4.465 

13.525 

94.425 

257.925 

501.175 

844.425 

1.320 

1.000 

13.375 

94.175 

256.725 

494.425 

842.075 

1,320 

0.825 

13.925 

94.425 

258.275 

508.625 

784.425 

1.320 

1.650 

14.425 

94.425 

256.575 

492.325 

782.925 

1.320 

2.475 

14.425 

91.675 

259.275 

488.425 

836.525 

1.320 

3.300 

15.275 

88.675 

262.675 

494.425 

770.625 

1.320 

4.125 

15.625 

87.625 

260.275 

494.425 

783.025 

1.320 

4.950 

15.825 

88.925 

252.525 

484.425 

835.225 

1.320 

5.775 

15.925 

92.075 

244.425 

483.675 

768.875 

1.320 

6.600 

16.025 

94.425 

244.425 

490.025 

765.025 

1.320 

7.425 

15.425 

87.675 

260.925 

500.625 

751.025 

1.360 

4.445 

13.425 

94.425 

256.425 

494.425 

844.425 

1.360 

1.000 
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75 


Mode  1 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Size 

Location 

11.925 

74.425 

222.325 

437.725 

657.225 

2.840 

0.825 

12.725 

77.625 

241.175 

379.375 

659.425 

2.840 

1.650 

13.425 

79.275 

218.825 

381.375 

661.075 

2.840 

2.475 

14.125 

79.425 

193.525 

432.425 

642.175 

2.840 

3.300 

14.425 

79.425 

186.825 

428.425 

642.925 

2.840 

4.125 

15.225 

79.425 

194.425 

370.625 

656.575 

2.840 

4.950 

15.575 

81.275 

221.075 

338.725 

637.625 

2.840 

5.775 

14.425 

79.425 

186.325 

439.425 

635.775 

2.880 

3.685 

12.125 

74.425 

228.025 

425.825 

657.125 

2.880 

1.000 

11.925 

73.875 

222.225 

433.525 

656.275 

2.880 

0.825 

12.725 

77.125 

240.525 

376.875 

657.125 

2.880 

1.650 

13.425 

78.925 

217.025 

381.275 

654.425 

2.880 

2.475 

14.125 

79.425 

191.675 

430.875 

640.025 

2.880 

3.300 

14.425 

79.425 

185.225 

424.425 

641.325 

2.880 

4.125 

15.225 

79.425 

194.425 

367.925 

650.075 

2.880 

4.950 

15.575 

81.075 

220.775 

336.925 

632.725 

2.880 

5.775 

14.375 

79.375 

184.425 

436.625 

635.775 

2.920 

3.665 

12.075 

74.125 

227.875 

421.425 

657.025 

2.920 

1.000 

11.925 

73.275 

222.125 

429.275 

655.575 

2.920 

0.825 

12.675 

76.625 

239.425 

374.425 

654.025 

2.920 

1.650 

13.425 

78.575 

215.175 

381.225 

644.425 

2.920 

2.475 

14.125 

79.325 

189.425 

429.175 

638.275 

2.920 

3.300 

14.425 

79.375 

183.675 

421.175 

639.425 

2.920 

4.125 

15.225 

79.425 

194.425 

365.325 

643.625 

2.920 

4.950 

15.575 

80.825 

220.425 

335.375 

627.925 

2.920 

5.775 

14.375 

79.175 

183.225 

433.425 

635.775 

2.960 

3.645 

12.075 

73.575 

227.775 

417.175 

657.025 

2.960 

1.000 

11.875 

72.725 

222.025 

425.025 

655.125 

2.960 

0.825 

12.675 

76.125 

238.925 

372.425 

650.625 

2.960 

1.650 

13.375 

78.225 

213.325 

381.175 

642.625 

2.960 

2.475 

14.075 

79.075 

187.925 

427.275 

636.775 

2.960 

3.300 

14.425 

79.225 

182.175 

417.575 

637.325 

2.960 

4.125 

15.175 

79.375 

194.425 

362.725 

637.225 

2.960 

4.950 

15.575 

80.625 

220.175 

333.825 

623.325 

2.960 

5.775 

14.325 

79.025 

181.725 

430.225 

635.775 

3.000 

3.625 

12.025 

73.025 

227.675 

412.925 

657.025 

3.000 

1.000 

11.875 

72.175 

221.925 

420.775 

654.425 

3.000 

0.825 

12.625 

75.625 

238.125 

370.525 

644.425 

3.000 

1.650 

13.375 

77.875 

211.425 

381.125 

636.875 

3.000 

2.475 

14.075 

78.875 

186.225 

425.175 

635.575 

3.000 

3.300 

14.425 

79.075 

180.675 

413.925 

634.425 

3.000 

4.125 

15.175 

79.225 

193.925 

360.225 

630.925 

3.000 

4.950 

15.575 

80.375 

219.425 

332.425 

618.875 

3.000 

5.775 

14.275 

78.825 

180.225 

426.925 

635.775 

3.040 

3.605 

12.025 

72.425 

227.525 

408.825 

656.925 

3.040 

1.000 

11.825 

71.625 

221.875 

416.525 

654.425 

3.040 

0.825 

12.625 

75.175 

237.175 

368.775 

643.075 

3.040 

1.650 

13.375 

77.525 

209.425 

381.075 

631.375 

3.040 

2.475 

14.025 

78.675 

184.425 

422.925 

634.425 

3.040 

3.300 

14.425 

78.925 

179.225 

410.325 

632.425 

3.040 

4.125 

15.175 

79.025 

193.075 

357.775 

624.425 

3.040 

4.950 

15.575 

80.175 

219.425 

331.125 

614.425 

3.040 

5.775 

14.275 

78.675 

178.775 

423.675 

635.725 

3.080 

3.585 

11.925 

71.925 

227.375 

404.425 

656.875 

3.080 

1.000 

11.825 

71.075 

221.775 

412.325 

654.425 

3.080 

0.825 

12.575 

74.425 

236.125 

367.175 

638.925 

3.080 

1.650 

13.325 

77.175 

207.775 

381.025 

626.175 

3.080 

2.475 

14.025 

78.425 

182.725 

420.425 

634.025 

3.080 

3.300 

14.425 

78.775 

177.825 

406.725 

629.425 

3.080 

4.125 

15.175 

78.875 

192.275 

355.425 

618.775 

3.080 

4.950 

15.525 

79.425 

219.325 

329.425 

610.525 

3.080 

5.775 

14.225 

78.425 

177.375 

420.325 

635.725 

3.120 

3.565 

11.925 

71.425 

227.175 

400.925 

656.675 

3.120 

1.000 
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Mode  1 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Size 

Location 

14.275 

75.725 

156.925 

335.825 

551.125 

4.000 

4.125 

14.425 

74.425 

183.425 

317.875 

524.425 

4.000 

4.950 

13.375 

74.275 

153.725 

344.425 

609.425 

4.040 

3.105 

11.375 

61.625 

208.075 

355.575 

593.175 

4.040 

1.000 

11.225 

60.525 

207.075 

352.675 

605.675 

4.040 

0.825 

12.025 

65.775 

194.425 

357.825 

544.425 

4.040 

1.650 

12.825 

71.075 

168.925 

344.425 

583.425 

4.040 

2.475 

13.575 

74.425 

151.425 

343.275 

605.575 

4.040 

3.300 

14.275 

75.575 

156.425 

333.425 

544.425 

4.040 

4.125 

14.425 

74.175 

183.275 

317.075 

522.725 

4.040 

4.950 

13.375 

74.075 

153.025 

342.025 

606.925 

4.080 

3.085 

11.375 

61.275 

206.575 

355.275 

589.025 

4.080 

1.000 

11.175 

60.175 

205.775 

352.125 

601.325 

4.080 

0.825 

11.925 

65.525 

194.425 

357.425 

544.425 

4.080 

1.650 

12.775 

70.875 

167.575 

344.425 

583.125 

4.080 

2.475 

13.575 

74.425 

150.575 

340.275 

602.225 

4.080 

3.300 

14.275 

75.425 

156.075 

331.225 

544.425 

4.080 

4.125 

14.425 

73.925 

183.175 

316.375 

520.875 

4.080 

4.950 

13.325 

73.925 

152.325 

339.225 

604.125 

4.120 

3.065 

11.325 

60.925 

205.025 

355.075 

584.425 

4.120 

1.000 

11.175 

59.425 

204.425 

351.625 

594.425 

4.120 

0.825 

11.925 

65.275 

194.425 

356.925 

544.425 

4.120 

1.650 

12.775 

70.725 

166.275 

344.425 

582.725 

4.120 

2.475 

13.525 

74.425 

144.425 

337.325 

594.425 

4.120 

3.300 

14.225 

75.275 

155.675 

329.025 

542.125 

4.120 

4.125 

14.425 

73.725 

183.025 

315.675 

519.175 

4.120 

4.950 

13.275 

73.725 

151.675 

336.425 

601.125 

4.160 

3.045 

11.325 

60.675 

203.525 

354.425 

581.025 

4.160 

1.000 

11.175 

59.425 

203.075 

351.225 

592.625 

4.160 

0.825 

11.925 

65.025 

192.825 

356.375 

544.425 

4.160 

1.650 

12.725 

70.575 

165.025 

342.525 

582.125 

4.160 

2.475 

13.525 

74.425 

144.425 

334.425 

594.425 

4.160 

3.300 

14.225 

75.125 

155.275 

326.925 

539.275 

4.160 

4.125 

14.425 

73.525 

182.925 

315.125 

517.675 

4.160 

4.950 

13.225 

73.525 

151.025 

333.775 

594.425 

4.200 

3.025 

11.275 

60.375 

201.925 

354.425 

577.125 

4.200 

1.000 

11.125 

59.225 

201.625 

350.925 

588.325 

4.200 

0.825 

11.925 

64.425 

191.125 

355.725 

543.675 

4.200 

1.650 

12.725 

70.425 

163.775 

340.075 

581.325 

4.200 

2.475 

13.425 

74.425 

144.425 

331.575 

591.575 

4.200 

3.300 

14.225 

74.425 

154.425 

324.425 

536.425 

4.200 

4.125 

14.425 

73.325 

182.775 

314.425 

516.325 

4.200 

4.950 

13.225 

73.325 

150.375 

331.125 

594.425 

4.240 

3,005 

11.275 

60.075 

200.325 

354.425 

573.375 

4.240 

1.000 

11.125 

58.925 

200.175 

350.725 

584.125 

4.240 

0.825 

11.925 

64.425 

189.425 

354,425 

542.825 

4.240 

1.650 

12.725 

70.275 

162.575 

337.575 

580.325 

4.240 

2.475 

13.425 

74.425 

144.425 

328.775 

587.825 

4.240 

3.300 

14.175 

74.425 

154.425 

322.925 

533.775 

4.240 

4.125 

14.425 

73.075 

182.625 

314.125 

515.075 

4.240 

4.950 

13.175 

73.125 

144.425 

328.575 

591.425 

4.280 

2.985 

11.275 

59.425 

194.425 

354.425 

569.425 

4.280 

1.000 

11.075 

58.625 

194.425 

350.525 

579.425 

4.280 

0.825 

11.875 

64.275 

187.775 

354.075 

542.125 

4.280 

1.650 

12.675 

70.125 

161.375 

335.125 

579.075 

4.280 

2.475 

13.425 

74.375 

144.425 

326.075 

584.025 

4.280 

3.300 

14.175 

74.425 

154.325 

321.075 

531.125 

4.280 

4.125 

14.425 

72.875 

182.425 

313.675 

514.025 

4.280 

4.950 

13.125 

72.925 

144.425 

326.075 

588.025 

4.320 

2.965 

11.225 

59.425 

194.425 

354.425 

566.275 

4.320 

1.000 

11.075 

58.325 

194.425 

350.425 

575.825 

4.320 

0.825 

11.875 

64.025 

186.125 

353.125 

541.575 

4.320 

1.650 

12.675 

69.425 

160.175 

332.675 

577.575 

4.320 

2.475 
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Mode  1 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Size 

Location 

13.425 

74.275 

144.425 

323.425 

580.175 

4.320 

3.300 

14.175 

74.425 

154.025 

319.275 

528.525 

4.320 

4.125 

13.075 

72.725 

144.425 

323.575 

584.425 

4.360 

2.945 

11.225 

59.225 

194.425 

354.425 

562.925 

4.360 

1.000 

11.025 

58.025 

194.425 

350.375 

571.875 

4.360 

0.825 

11.825 

63.825 

184.425 

352.075 

541.125 

4.360 

1.650 

12.625 

69.425 

159.075 

330.175 

575.825 

4.360 

2.475 

13.425 

74.175 

144.425 

320.825 

576.275 

4.360 

3.300 

14.125 

74.325 

153.775 

317.575 

526.075 

4.360 

4.125 

13.025 

72.525 

144.425 

321.175 

580.925 

4.400 

2.925 

11.175 

58.925 

193.875 

354.425 

559.425 

4.400 

1.000 

11.025 

57.775 

194.125 

350.325 

568.025 

4.400 

0.825 

11.825 

63.625 

182.875 

350.925 

540.825 

4.400 

1.650 

12.625 

69.425 

157.925 

327.775 

573.875 

4.400 

2.475 

13.375 

74.075 

144.425 

318.325 

572.325 

4.400 

3.300 

14.125 

74.175 

153.575 

315.925 

523.625 

4.400 

4.125 

13.025 

72.325 

144.425 

318.825 

577.225 

4.440 

2.905 

11.175 

58.675 

192.225 

354.425 

556.775 

4.440 

1.000 

10.925 

57.425 

192.575 

350.325 

564.275 

4.440 

0.825 

11.775 

63.425 

181.275 

344.425 

540.675 

4.440 

1.650 

12.575 

69.425 

156.825 

325.325 

571.625 

4.440 

2.475 

13.375 

73.925 

143.775 

315.875 

568.375 

4.440 

3.300 

14.125 

73.925 

153.375 

314.425 

521.325 

4.440 

4.125 

12.925 

72.125 

144.425 

316.525 

573.425 

4.480 

2.885 

11.125 

58.425 

190.575 

354.425 

553.925 

4.480 

1.000 

10.925 

57.225 

191.025 

350.325 

560.725 

4.480 

0.825 

11.775 

63.225 

179.425 

344.425 

540.525 

4.480 

1.650 

12.575 

69.425 

155.775 

322.925 

569.175 

4.480 

2.475 

13.375 

73.875 

143.175 

313.425 

564.425 

4.480 

3.300 

14.075 

73.825 

153.175 

312.925 

519.075 

4.480 

4.125 

12.925 

71.925 

144.425 

314.275 

569.425 

4.520 

2.865 

11.125 

58.175 

188.925 

354.425 

551.325 

4.520 

1.000 

10.925 

56.925 

189.425 

350.325 

557.325 

4.520 

0.825 

11.725 

63.025 

178.175 

344.425 

540.425 

4.520 

1.650 

12.575 

69.325 

154.425 

320.525 

566.525 

4.520 

2.475 

13.325 

73.825 

142.575 

311.125 

560.425 

4.520 

3.300 

14.075 

73.675 

153.025 

311.575 

516.925 

4.520 

4.125 

12.875 

71.725 

144.425 

312.075 

565.875 

4.560 

2.845 

11.075 

57.925 

187.325 

354.325 

544.425 

4.560 

1.000 

10.925 

56.725 

187.925 

350.275 

554.075 

4.560 

0.825 

11.725 

62.825 

176.675 

344.425 

540.375 

4.560 

1.650 

12.525 

69.225 

153.725 

318.175 

563.675 

4.560 

2.475 

13.325 

73.725 

142.025 

308.875 

556.525 

4.560 

3.300 

14.075 

73.425 

152.875 

310.275 

514.425 

4.560 

4.125 

12.825 

71.425 

144.425 

309.425 

561.925 

4.600 

2.825 

11.075 

57.725 

185.725 

354.125 

544.425 

4.600 

1.000 

10.875 

56.425 

186.325 

350.225 

551.025 

4.600 

0.825 

11.725 

62.625 

175.175 

343.925 

540.325 

4.600 

1.650 

12.525 

69.125 

152.725 

315.825 

560.625 

4.600 

2.475 

13.275 

73.625 

141.425 

306.675 

552.575 

4.600 

3.300 

14.025 

73.275 

152.725 

309.025 

512.875 

4.600 

4.125 

12.825 

71.275 

144.425 

307.875 

558.075 

4.640 

2.805 

11.025 

57.525 

184.125 

353.775 

544.425 

4.640 

1.000 

10.875 

56.225 

184.425 

350.125 

544.425 

4.640 

0.825 

11.675 

62.425 

173.675 

342.325 

540.225 

4.640 

1.650 

12.425 

69.025 

151.775 

313.525 

557.375 

4.640 

2.475 

13.275 

73.525 

140.925 

304.425 

544.425 

4.640 

3.300 

14.025 

73.125 

152.625 

307.925 

511.025 

4.640 

4.125 

12.775 

71.075 

144.425 

305.875 

554.125 

4.680 

2.785 

11.025 

57.275 

182.525 

353.375 

542.875 

4.680 

1.000 

10.875 

56.025 

183.225 

344.425 

544.425 

4.680 

0.825 

11.675 

62.275 

172.225 

340.625 

540.025 

4.680 

1.650 

12.425 

68.925 

150.825 

311.275 

554.025 

4.680 

2.475 

Mode  1 

13.275 
14.025 

12.725 

10.925 

10.825 

11.625 

12.425 

13.225 

13.925 

12.675 

10.925 

10.825 

11.625 

12.425 

13.225 

13.925 

12.625 

10.925 

10.775 

11.575 

12.375 

13.175 

13.925 

12.625 

10.925 

10.775 

11.575 

12.375 

13.175 

13.925 

12.575 

10.925 

10.725 

11.525 

12.375 

13.175 

13.925 

12.525 

10.875 

10.725 

11.525 

12.325 

13.125 

13.925 

12.425 

10.875 

10.675 

11.425 

12.325 

13.125 

13.875 

12.425 

10.825 

10.675 

11.425 

12.275 
13.075 

13.875 

12.375 

10.825 

10.625 

11.425 

12.275 


Mode  2 

73.425 

72.925 

70.875 
57.075 

55.775 

62.125 

68.825 

73.325 

72.775 

70.675 

56.875 

55.575 

61.925 

68.725 

73.225 

72.575 

70.425 

56.675 

55.375 

61.825 

68.625 

73.125 

72.375 

70.225 

56.425 

55.175 

61.675 

68.575 
73.025 

72.175 
70.025 

56.275 

54.425 

61.525 

68.425 

72.925 
72.025 

69.425 

56.125 

54.425 

61.375 

68.375 

72.825 

71.825 

69.425 

55.925 

54.425 

61.275 

68.325 

72.725 

71.625 

69.375 

55.775 

54.425 

61.125 

68.225 

72.625 

71.425 

69.175 

55.575 

54.225 
61.025 

68.175 


Mode  3 

140.525 

152.525 

144.275 

180.925 

181.625 

170.825 

144.425 
140.075 

152.425 

143.825 

179.325 
180.075 

169.425 

144.425 

139.425 

152.325 

143.425 

177.775 

178.575 
168.075 

144.425 

139.225 

152.275 
143.025 

176.225 
177.025 

166.725 

144.425 

138.825 

152.225 

142.675 

174.425 

175.525 

165.375 

144.425 

138.425 

152.125 

142.275 

173.225 
174.025 
164.075 

144.425 

138.125 
152.075 

141.925 

171.725 

172.525 

162.825 

144.425 

137.825 
152.025 

141.575 

170.275 
171.075 

161.575 

144.125 

137.525 

151.925 

141.275 

168.825 

169.425 

160.375 

143.425 


Mode  4 

302.525 

306.875 

303.875 

352.925 

344.425 

338.925 
309.025 

300.525 

305.875 

301.925 

352.325 

344.425 

337.125 

306.825 

294.425 

304.425 

300.125 

351.675 

344.425 

335.325 

304.425 

294.425 

304.175 

294.425 

350.875 

344.425 

333.425 

302.525 

294.425 

303.425 

294.425 

344.425 

344.425 

331.575 

300.425 

293.275 

302.725 

294.425 

344.425 

344.425 

329.425 

294.425 

291.625 

302.125 

293.225 

344.425 

344.425 

327.775 

294.425 
290.075 

301.625 

291.625 

344.425 

344.425 

325.825 

294.425 

288.525 

301.125 
290.025 

344.425 

344.425 

323.875 

292.525 


Mode  5 

544.425 

509.225 

550.175 

541.325 
543.075 

539.425 

550.525 

540.925 

507.575 

544.425 

539.425 

540.825 

539.425 

544.425 
537.075 

505.925 

542.225 

538.725 

538.775 

538.925 

543.175 

533.325 

504.425 

538.225 

537.725 

536.925 

538.225 

539.425 

529.425 

503.125 

534.275 

536.875 

535.375 

537.375 

535.575 

525.925 

501.825 

530.325 

536.225 

533.925 

536.375 

531.675 

522.325 

500.675 

526.375 

535.675 

532.725 

535.175 

527.725 

518.825 

494.425 

522.425 

535.325 

531.725 

533.725 

523.775 

515.325 

494.425 

518.575 
535.025 

530.925 

532.125 

519.425 


Size 

4.680 

4.680 

4.720 

4.720 

4.720 

4.720 

4.720 

4.720 

4.720 

4.760 

4.760 

4.760 

4.760 

4.760 

4.760 

4.760 

4.800 

4.800 

4.800 

4.800 

4.800 

4.800 

4.800 

4.840 

4.840 

4.840 

4.840 

4.840 

4.840 

4.840 

4.880 

4.880 

4.880 

4.880 

4.880 

4.880 

4.880 

4.920 

4.920 

4.920 

4.920 

4.920 

4.920 

4.920 

4.960 

4.960 

4.960 

4.960 

4.960 

4.960 

4.960 

5.000 

5.000 

5.000 

5.000 

5.000 

5.000 

5.000 

5.040 

5.040 

5.040 

5.040 

5.040 


Location 

3.300 

4.125 

2.765 

1.000 

0.825 

1.650 

2.475 

3.300 

4.125 

2.745 

1.000 

0.825 

1.650 

2.475 

3.300 

4.125 

2.725 

1.000 

0.825 

1.650 

2.475 

3.300 

4.125 

2.705 

1.000 

0.825 

1.650 

2.475 

3.300 

4.125 

2.685 

1.000 

0.825 

1.650 

2.475 

3.300 

4.125 

2.665 

1.000 

0.825 

1.650 

2.475 

3.300 

4.125 

2.645 

1.000 

0.825 

1.650 

2.475 

3.300 

4.125 

2.625 

1.000 

0.825 

1.650 

2.475 

3.300 

4.125 

2.605 

1.000 

0.825 

1.650 

2.475 
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Mode  1 

10.175 

10.925 

11.825 

12.675 

11.425 

10.325 

10.175 

10.925 

11.775 

12.675 

11.425 

10.325 

10.125 

10.925 

11.775 

11.375 

10.275 

10.125 

10.925 

11.725 

11.325 

10.275 
10.075 

10.875 

11.725 

11.325 

10.225 
10.075 

10.875 

11.725 

11.275 

10.225 
10.075 

10.825 

11.675 

11.225 

10.175 
10.025 

10.825 

11.675 

11.175 

10.175 
10.025 

10.825 

11.625 

11.125 

10.125 

9.425 

10.775 

11.625 
11.075 

10.125 

9.425 

10.775 

11.625 

11.075 

10.125 

9.425 

10.725 

11.575 
11.025 
10.075 

9.425 


Mode  2 

51.725 

59.375 

66.775 

69.425 

64.225 

53.175 

51.675 

59.325 

66.725 

69.375 

63.925 

53.125 

51.575 

59.275 

66.625 

63.775 
53.075 

51.525 

59.225 

66.575 

63.525 
53.025 

51.425 

59.175 

66.425 

63.275 

52.925 

51.375 

59.175 

66.375 
63.075 

52.875 

51.325 

59.125 

66.325 

62.825 

52.825 

51.275 
59.075 

66.225 

62.575 

52.775 

51.225 
59.025 

66.125 

62.375 

52.725 

51.175 

58.925 
66.025 

62.125 

52.725 

51.125 

58.925 

65.925 

61.875 

52.675 
51.075 

58.925 

65.825 

61.675 

52.625 
51.025 


Mode  3 

143.925 

140.575 

133.425 
135.025 

135.425 

142.925 
143.025 

139.425 

133.175 
135.025 

135.275 
142.025 
142.075 

139.275 

132.925 
135.025 

141.125 

141.125 

138.625 

132.775 

134.425 

140.275 

140.225 
138.025 

132.575 

134.425 

139.425 

139.375 

137.425 

132.375 

134.425 

138.625 

138.425 

136.875 

132.225 

134.175 

137.825 

137.675 

136.375 
132.075 

133.925 
137.075 

136.825 

135.875 

131.925 

133.775 

136.325 
136.075 

135.375 

131.825 

133.575 

135.625 

135.275 

134.425 

131.725 

133.375 

134.425 

134.425 

134.425 

131.625 

133.175 

134.225 

133.825 


Mode  4 
308.075 

284.425 

262.825 

271.525 

266.775 

303.425 

305.925 

283.225 

261.925 

271.375 

266.125 

301.325 

303.825 

281.675 
261.075 

265.425 

294.425 

301.725 

280.125 

260.275 

264.425 

294.425 

294.425 

278.625 

259.425 

264.275 

294.425 

294.425 

277.175 

258.825 

263.725 

293.225 

294.425 

275.725 

258.175 

263.175 

291.225 

293.425 

274.325 

257.575 

262.675 

289.275 

291.425 

272.925 

256.925 

262.175 

287.375 

289.425 

271.625 

256.425 

261.675 

285.425 

287.425 

270.325 
256.025 

261.175 

283.575 

285.425 
269.075 

255.575 

260.725 

281.775 

283.525 


Mode  5 

519.125 

471.225 

444.425 

464.125 

444.425 

513.425 

517.175 

467.925 

442.275 

463.125 

442.925 

510.725 
515.075 

464.425 

439.425 

440.425 

507.925 

512.825 

461.525 

436.875 

437.925 
505.025 

510.425 

458.325 

434.325 

435.575 
502.025 

507.925 

455.175 

431.875 

433.225 

494.425 

505.225 
452.075 

429.425 

430.925 

494.425 

502.425 

444.425 

427.225 

428.775 

492.725 

494.425 

444.425 
425.025 

426.625 

489.425 

494.425 

442.925 

422.925 

424.425 

486.325 

493.625 
440.075 

420.925 

422.575 

483.125 

490.575 

437.175 
419.075 

420.675 

479.425 

487.425 


Size 

5.880 

5.880 

5.880 

5.880 

5.920 

5.920 

5.920 

5.920 

5.920 

5.920 

5.960 

5.960 

5.960 

5.960 

5.960 

6.000 

6.000 

6.000 

6.000 

6.000 

6.040 

6.040 

6.040 

6.040 

6.040 

6.080 

6.080 

6.080 

6.080 

6.080 

6.120 

6.120 

6.120 

6.120 

6.120 

6.160 

6.160 

6.160 

6.160 

6.160 

6.200 

6.200 

6.200 

6.200 

6.200 

6.240 

6.240 

6.240 

6.240 

6.240 

6.280 

6.280 

6.280 

6.280 

6.280 

6.320 

6.320 

6.320 

6.320 

6.320 

6.360 

6.360 

6.360 


Location 

0.825 

1.650 

2.475 

3.300 

2.165 

1.000 

0.825 

1.650 

2.475 

3.300 

2.145 

1.000 

0.825 

1.650 

2.475 

2.125 

1.000 

0.825 

1.650 

2.475 

2.105 

1.000 

0.825 

1.650 

2.475 

2.085 

1.000 

0.825 

1.650 

2.475 

2.065 

1.000 

0.825 

1.650 

2.475 

2.045 

1.000 

0.825 

1.650 

2.475 

2.025 

1.000 

0.825 

1.650 

2.475 

2.005 

1.000 

0.825 

1.650 

2.475 

1.985 

1.000 

0.825 

1.650 

2.475 

1.965 

1.000 

0.825 

1.650 

2.475 

1.945 

1.000 

0.825 


APPENDIX  C 

INPUT  AND  OUTPUT  FOR  THE  FINITE  ELEMENT  MODEL 
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The  following  appendix  contains  the  inputs  and  outputs  of  two  separate  runs  of  a 
finite  element  model  discussed  earlier  in  the  text  in  reference  [32].  The  finite  element 
code  uses  English  units  so  the  inputs  were  converted  from  the  MKS  units  used  in  the  text 
to  values  listed  below.  The  finite  element  code  also  calculated  72  modal  frequencies  but 
only  the  first  six  frequencies  are  listed  in  this  appendix.  Both  runs  were  for  glass-epoxy 
cantilever  composite  beams  with  identical  lamina  properties;  however,  the  beams  had 
different  dimensions.  In  the  first  run,  the  beam  was  2.54  cm  (1”)  wide  and  0.18  cm 
(0.072”)  thick.  The  second  beam  was  2.33  cm  (0.917”)  wide  and  0.15  cm  (0.061”)  thick. 


Variable  definition 


NPLY  =  number  of  lamina 
NEM  =  number  of  elements  in  the  mesh 
NDF  =  number  of  degrees  of  freedom  per  node 
TH(I)  =  lamina  orientation  (degrees) 

AL  =  beam  length  (inches) 

XO  =  location  of  beam  clamp  (inches) 

HT  =  beam  thickness  (inches) 

B  =  beam  width  (inches) 

V12  =  Poisson  ration 

El,  E2  =  principal  moduli  (psi) 

G12,  G13,  G23  =  shear  moduli  (psi) 

RHO  =  density  (lb-secVin‘‘) 
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Run  1 
Inputs: 

NPLY  =  8 
NEM=  18 
NDF  =  4 
AL=  10.25 
XO  =  0.0 
HT  =  0.072 
B  =  1.0 
V12  =  0.27 
El  =6.14E+06 
E2=  1.7E+06 
G12=  1.03E+06 
G13  =  1.03E+06 
G23  =  0.93E+06 
RHO  =  1.805E-04 

TH(I)  0.0  90.0  0.0  90.0  90.0  0.0  90.0  0.0 

Outputs; 

Mode  1  Mode  2  Mode  3  Mode  4  Mode  5 

11.29  Hz  71.01Hz  199.14Hz  394.93  Hz  659.76  Hz 


Run  2 
Inputs: 

NPLY  =  8 
NEM=  18 
NDF  =  4 
AL  =  10.25 
XO  =  0.0 
HT  =  0.061 
B  =  0.917 
V12  =  0.27 
El  =  6.14E+06 
E2=  1.7E+06 
G12=  1.03E+06 
G13  =  1.03E+06 
G23  =  0.93E+06 
RHO  =  1.805E-04 

TH(I)  0.0  90.0  0.0  90.0  90.0  0.0  90.0  0.0 

Outputs: 

Mode  1  Mode  2  Mode  3  Mode  4  Mode  5 

15.31Hz  96.23  Hz  270.87  Hz  534.90  Hz  893.14  Hz 


Mode  6 
998.54  Hz 


Mode  6 
1350.89  Hz 
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